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Preface 



One hundred years ago (1904) Hermann Minkowski [58] posed a problem: to re- 
construct an even function / on the sphere S 2 from knowledge of the integrals 

Mf(C)= [ fds 
Jc 

over big circles C. Paul Funk found an explicit reconstruction formula for / from 
data of big circle integrals. Johann Radon studied a similar problem for the Eu- 
clidean plane and space. The interest in reconstruction problems like Minkowski- 
Funk’s and Radon’s has grown tremendously in the last four decades, stimulated by 
the spectrum of new modalities of image reconstruction. These are X-ray, MRI, 
gamma and positron radiography, ultrasound, seismic tomography, electron mi- 
croscopy, synthetic radar imaging and others. The physical principles of these 
methods are very different, however their mathematical models and solution meth- 
ods have very much in common. The umbrella name reconstructive integral geom- 
etry 1 is used to specify the variety of these problems and methods. 

The objective of this book is to present in a uniform way the scope of well- 
known and recent results and methods in the reconstructive integral geometry. 
We do not touch here the problems arising in adaptation of analytic methods to 
numerical reconstruction algorithms. We refer to the books [61], [62] which are 
focused on these problems. 

Various aspects of interplay of integral geometry and differential equations 
are discussed in Chapters 7 and 8. The results presented here are partially new. 

The book is an extended version of a lecture course which was read for 
students of Tel Aviv University. Chapter 1 contains a mini-course in distribution 
theory, harmonic analysis and geometry. Not much of this knowledge is necessary 
for reading and understanding Chapters 2, 3 and 6. 



^laschke’s term “Integralgeometrie” looks somehow redundant; the word geometry itself has, 
since ancient Greek times, meant calculation of lengths, areas, volumes, i.e., some integrals. 




Notation 



Z - the algebra of integers, 

R - the field of real numbers, 

R+ - the set of non-negative numbers 
C - the field of complex numbers 
j = 2tti, i = yJ~—\ 

f+ = max (/, 0) , f—= max (-/, 0) 

r (A) , A G C - Gamma-function 

V - a vector space over M of finite dimension 

V* - the space dual to a vector space V 

dx or dV - the volume form in a coordinate space V 

E or E n - a Euclidean space of dimension n < oc 

D (V) - the space of smooth functions of compact support (test functions) in the 
vector space V 

K (V) - the space of smooth densities of compact support (test densities) in V 
S (V) - the space of smooth fast decreasing functions in V 

F (/) = / - the Fourier transform of a function / or of the density / dx defined 
in V; / (£) = / v exp (-j£r) / (x) dz 

F* (/) = / v * exp (j£x) / (£) d£ - the adjoint Fourier transform 




Chapter 1 

Distributions and Fourier Transform 



1.1 Introduction 

The Riemann theory of integration in a real vector space of finite dimension is 
sufficient for our purposes except for Section 1.7, where we outline the Fourier- 
Plancherel theory. Let V be a real vector space of finite dimension n. Choose a 
coordinate system x = (aq, . . . , x n ) and take the measure dx = dxi . . . dx n in V. 
A function / : V — > C is called Riemann integrable , if it is continuous almost 
everywhere in V, (i.e., except for a zero measure set) and the norm 

ll/lli = sup [ \ft\dx 

B,t JB 

is finite, where B denotes an arbitrary ball in V and f t = f if |/| < t and f t = 0 
otherwise (truncation of /). We set 




where B ( r ) is the ball of radius r centred at the origin (cubes can be taken instead 
of balls etc.). We shall use standard tools of Riemann integration theory including 
Fubini’s Theorem. The full range Fubini theorem will be stated in Section 1.7 in 
the framework of Lebesgue integration theory. 

1.2 Distributions and generalized functions 

Let V be a finite-dimensional vector space, U an open set in V. A function / : 
U — > C is called smooth or of class (7°°, if it has all derivatives with respect to 
a system of coordinates in V. This property is equivalent to the condition: / has 
a differential of arbitrary order at each point x € V and is preserved under any 
linear (or smooth) coordinate change. 

Definition. The notation D = D (U) means the space of smooth functions (j) in V 
with compact support supp (j) (e U (called test functions) . This is a vector space 
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over the field C of complex numbers. The natural convergence in D ( U ) is defined 
as follows: {fa, k = 1,2,...} converges to </>, if 

(i) U supp fa C K <^U and 

(ii) D q (fa — (/))=> 0 uniformly in U for arbitrary q = where 

^lH 1 -Qn 

D q = 

(dxi)* 1 . . . (dx n Y n 

and x\ , . . . , x n is a linear coordinate system in V. A functional u : D (U) — > C 
is called continuous with respect to the natural convergence, if u (fa) — > u (</>) 
for any convergent sequence fa — > 0. Any linear continuous functional on 
D (V) is called Schwartz’s distribution in U. The space of all distributions is 
denoted D' ( U ). 

A test density in U is a smooth density p with compact support. For a choice 
of the coordinate system we can write p = <fi dx where dx = dxi A • • • A dx n and 
(j) = pj dx. The space of test densities with support in U is denoted K ( U ). This 
space is supplied with the natural convergence: a sequence pk converges to p as the 
sequence of test functions fa — pk/dx converges to (f) = p/dx in the sense (i)— (ii). 
A linear continuous functional on the space K (U) is called Sobolev’s generalized 
function in U. 1 The linear operations, multiplication by a smooth function and 
derivatives are well defined in the space D' ( U ) of distributions and in the space 
K' ( U ) of generalized functions. 

Example 1. A function / : V — » C is called locally integrable, if its restriction to 
an arbitrary compact set K C U is integrable. Any locally integrable function / 
defines a generalized function [/] as 

[/] (p) = f fp = j p = (jxlx eK(U). 

The functional [/] equals zero if and only if the function / vanishes almost every- 
where. The density fdx defines a distribution denoted [/dx]. 

Example 2. A Dirac function at a point aG Vis the functional 5 a (4> dx) = (p(a) on 
the space K (V). Dirac-distribution S a dx is the continuous functional on D (V): 
8 a dx (fa = (j) (a). 

Example 3. For a function </> : R — > C the limit 







(j)dx 

x 



is called the principal value of the divergent integral. The limit exists for any 
integrable function (j) with compact support that has derivative at the point x = 0. 
The functional [dx/x] is linear and continuous in D (M), i.e., a distribution in R. 

x de Rham’s construction of currents [18] contains generalized functions (currents of degree 
0) and distributions (currents of degree n) as particular cases. 
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Example 4. Euler kernels. For an arbitrary complex A, Re A > 0, we define the 
functional on S (R), 

flA w = iMr <*> = m l xA_V(x,dI ' 

where x A_1 = exp ((A — 1) lnx), lnx G R for x > 0. This is a tempered distribution 
which depends analytically on A, i.e., H x (</>) is a holomorphic function of A for 
each (j). 

Proposition 1.1 The family H x has holomorphic continuation on the whole com- 
plex plane C. We have H~ k = <S<*> (0) for k = 0, 1, 2, . . . 

◄ We have for Re A > 1, for </>' = d</>/dx, 

(H x (VO)' = -H x (</>') = J x^ 1 ft (x)dx = r (fz T) j xX ~ 2< f > ( x ) dx ^ 

where we integrate by parts and take into account that x x ~ 1 <f(x) vanishes at 
x = 0 and at x = oo. The right side is equal to H x ~ l (4>) and the differential 
equation ( H x )' = # A_1 holds. Now we can define H x = (i/ A+1 ) / for Re A > — 1. 
This formula defines analytic family H x which coincides with the Euler family for 
Re A > 0. Next we can extend this family for Re A > — 2 and so on. The union 
of these continuations gives the holomorphic family in the whole complex plane. 
This continuation fulfils the differential equation. Calculate H°; take a function 
'ipo £ S (R) that coincides with exp (— x) for x > 0. We have for Re A > 0 

H x {ft) = H x (<t>-4> (0) VO + V> (0) H x (VO , 

H x (VO = f x x_1 exp (-x) dx = 1, 

H x (4>-<t> (0) VO = f (<f> ( X ) - V> (0) v> (x)) dx. 

The last integral has analytic continuation for Re A > — 1 since the function (x) — 
4> (0) (x) vanishes at x = 0. On the other hand the dominator T has a pole at 
A = 0 and the right sides vanishes at A = 0. Therefore H x (</>) — » 0 (0) as A — > 0, 

i.e., H° = (So- From the differential equation we find H~ k = (i/ 0 )^ = dQ k \ 
k = 1 , 2 ,.... ► 

We define H x = x A_1 dx/T (A) = (— x)+ _1 d x/T (A) for Re A > 0 and extend 
for all A by H x ( (j) ) = H x (if) where (x) = (j) (— x). 

Example 5. Consider the Cauchy integral 



lim^ I (x ± ei) x cp(x) dx 
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for a test density ipdx G S( R). These limits exists for any A G C and defines two 
families of generalized functions in R denoted by ( x =b 0z) A , respectively. Indeed, 
integrating by parts k times yields 

J (x±ei) x (p{x)dx = + ^(A + fc) f ( x±£t ) X+k V {k) (x) d®. 

The right-hand side has meromorphic continuation at the half-plane Re A > — k — 1 
with occasional poles at A = -1,-2,...,— k where the dominator vanishes. The 
numerator vanishes too which can be seen from the partial integration in the 
opposite direction. Therefore the right-hand side has holomorphic continuation to 
the half-plane and hence, to the whole plane. These generalized functions depend 
analytically on the parameter A G C since the functions (x ± ei) x do. We have 
d(x±0i) x /dx = A(x±0z) a_1 . The products (x±(h) A dx are tempered distributions. 
Two important equations hold: 



(x + (h) 1 dx + (x — Oz) 1 da; = 2 



dx 

x 



(x - Ot) 1 



(x + 0?) = j^o- (1-1) 



Moreover, for any A^— 1,-2,. ..we have 

( x ± 0z) A = x A -b exp (±\ni) x A . (1.2) 



This equation is obvious for the case Re A > -T where both sides are locally 
integrable. It holds for all A ^ —k due to uniqueness of analytic continuation. 
Operations. Linear operations: addition and multiplication by scalars are well 
defined in the space of distributions D' (U) and in the space of generalized func- 
tions K' (U) as operations on linear functionals. For any smooth function a in U 
the multiplication operator u i— > au is also well defined by au ((/)) — u (u</>), where 
(j) is any test function (density) in U. Partial derivatives are defined by the formula 



9U (A\ 




which is consistent with differentiation of smooth densities u = f dx. Let A = 
A (x, D) = Y^Ai (x) D l be an arbitrary linear differential expression with smooth 
coefficients Ai in U C V. It defines a linear operator on distributions u and on 
generalized functions v by 



Au {(j)) = u (A*<f )) , Av (</>dx) = v (A* (</>) dx) , 



where A* is the (formal) adjoint operator. The definition of the operator A* de- 
pends on the choice of volume density dx and can be read from the identity 



J A*(fyij ) dx = J 



<^A^dx, 
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where </>, i/j are arbitrary test functions in U. An independent definition of action 
of a vector field T to a generalized function v can be given as 

T (v) ( p ) = v (T* (p)) , T* ( p ) = -L ( T ) p, 

where L(T) p is action of the Lie derivative along T on the density p, see Sec- 
tion 1.6. 

Let A be again a linear differential operator in an open set U and s be a 
point in U. A fundamental solution for A with source at a point s is a generalized 
function F s (or the distribution F s dx ) in U that satisfies the equation AF S = 5 S 
where S s is the Dirac function at s. 

1.3 Tempered distributions 

Consider a vector space V of dimension n. The Schwartz space S (V) of test 
functions consists of smooth functions <p in V that satisfy the inequality 

\x p D q (p(x)\ <C(p,q ), x p = x Pl . . . x^ 71 (1.3) 

for any vectors p,q e Z n and some constant C (p, q) in a linear coordinate system 
x = (#i, . . . , x n ). The convergence in S (V) is defined by means of the system of 
norms 

IMIfc = . .max sup \x p D q ip(x) \ ,k = 0, 1,2, , 

M + M<fcx€V 

so that (fi —> (p if and only if || — ip\\ k — > 0 as i — > oo for any k = 0,1, — 
The convergence does not depend on the choice of a coordinate system. The space 
S(V) is the basis for L. Schwartz’s theory of the Fourier transform of distributions. 
The dual space S'(V) is the space of all linear continuous functions on ^(V); an 
element of this space is called a tempered distribution on V. 

A Schwartz density in V is by definition a product p = (f) dx where </> is 
a Schwartz function, x = (aq, . . . , x n ) is a linear coordinate system and dx = 
dx\ A • • • Adx n is a volume form in V 2 . We use the notation S(V) dx for the space 
of Schwartz densities in V. An arbitrary element of the dual space (S(V)dx) f is 
called a tempered generalized function on V. For an arbitrary tempered generalized 
function v the product u = vdx is a tempered distribution and any tempered 
distribution can be written in this form. 

Examples 6. If a function / in V has at most polynomial growth at infinity, i.e., 

| / (a) | < C (\x\ + l) q for some constants <7, C, then the generalized function [/] is 
tempered. The distributions considered in Examples 2-5 are tempered. 

Fourier transform. We fix a coordinate system x = (aq,...,x n ) in V; it 
generates the bijection V = M n and the measure (volume density) dx = dxi A 

2 It is easy to give a definition of Schwartz densities which does not depend on the choice of 
the specific volume form da:. 
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• • • A dx n . For any p > 1 the space L P (V) is constructed by means of the measure 
dx in V. The Fourier transform of an integrable function / is given by the integral 

/(&,•••,£«) = / /(aJi,...,a:n)exp(-j(^ia;i + ---+C n a: n ))dx. 

Jv 

We interpret the vector £ = (£i, . . . , £ n ) as a point in the dual space V* in such a 

way that the bilinear form (£,x) = £ 1 X 1 H + £ n x n is the value of a function £ 

on the vector x. We call the coordinate system £i, . . . ,£ n dual to the coordinate 
system xi, . . . , x n and write the Fourier integral in the abbreviated form 

F(fdx)(£) = /(£) = J /(x)exp(-j {£,x))dx. 



The improper Riemann integral exists for each point £ G V* and the function / 
is continuous and bounded: 

1/(01 < ll/lli- (1-4) 

We write the Fourier integral as an operator applied to the density / dx, its image 
is a function. This makes the operator F independent of the choice of coordinates 
in V. The adjoint Fourier integral is applied to a function and is equal to a density. 
Therefore we write it in the form 

F*(g)(x) = ^ J <7(0 ex P(j (00)d£j dx, 

where d£ = d£i A • • • Ad£ n for the dual coordinate system £i, . . . , £ n . The operators 
F and F* are inverse one to another; it will be stated in Section 1.7. Note that 
the operator F * also does not depend on the choice of linear coordinates x in V. 
Indeed, for another choice of coordinates y we have y = Ax where A is a non- 
singular matrix. The dual system coordinates p in V* relate to £ by the formula 
py = £x which implies p = A' -1 £, where A! is a transposed matrix. This yields 
dy = |det A\ dx and dp = |detx4| -1 d£. It follows that dpdy = d£dx and 



exp (j < p , y)) (p) $ (p) dp dy = 



exp (j (£, x)) (£) ^ (£) d£ dx, 



q.e.d. Nevertheless, to simplify our notation we shall frequently omit the factor 
dx. 

Basic properties. 1. If / is an integrable function and possesses integrable 
derivative df /dx* for some i, 1 < i < n, then 



F 



dxi 



j iihi). 



2. If / and x*/ are integrable functions, then the function F (/) has a continuous 
derivative in £* and 



d F(f) 

dti 



-j F(xif). 
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3. The Fourier transform acts in the Schwartz spaces. Recall that S (V) dx means 
the space of Schwartz test densities in V. 

Proposition 1.2 The Fourier transform acts continuously in the Schwartz spaces: 

F : S(V ) dx -> S(V*). 

◄ It is well defined, since S(V) C L\. Now we state the inclusion F(S (V) dx) C 
S (V*). Take an arbitrary function cp G S. Any derivative (p^ belongs to 5, hence 
by the property 1 the function (p is equal to 0(|£| -9 ) for arbitrary q. On the 
other hand the function x p (p(x) is again in S for any p. By 2, the function (p has 
derivatives of arbitrary order and any derivative is bounded in V*. Moreover, for 
arbitrary p,q we have 

(_l)|g|jkl-|p| j p(£ ) p- c ^) = ^D\F{(p). 

The left side is the Fourier image of a function which is a linear combination of 
functions x k (p^ l \ ki < qi,k < pi. These functions belong to Li, hence the left side is 
bounded by (1.4). The equation implies that the function F(<p) belongs to 5(V / ). It 
is easy to check by means of this equation that the operator F : 5(V)dx — » S(V*) 
is continuous. The same arguments are valid for the operator F*. ► 

Proposition 1.2 implies that the adjoint Fourier transform makes a test den- 
sity of a test function and defines a continuous operator F* : S(V*) — > *S'(V)dx. 

Definition. The Fourier transform of tempered distributions is the dual operator 

F f : S'(V) -> (S(V*)d€)' 

defined by F f (v)((p) = v(F(ip)). The function F' (v) is a tempered generalized func- 
tion. The adjoint Fourier transform is defined in a similar way F*' : (S(V)dx) / —> 
5'(V*); it makes a tempered distribution of a tempered generalized function. We 
call F f and F*' Fourier- Schwartz transforms. They are inverse one to another. 

Proposition 1.3 The Fourier- Schwartz transforms match the Fourier and adjoint 
Fourier transforms of test functions and make the commutative diagram 



F' : S'(V) 


- (S(V*)dO' 


F*' : (5(V*) d£)' - 


- S'(V) 


U 


u 


U 


u 


F: S (V) dx 


-> S (V*) 


F*: S (V*) 


S (V) dx. 



◄ Indeed, any test density g = f dx defines a tempered distribution by means 
of integration [g] (cp) = J f(pdx. Then F' (g) is the tempered generalized function 
in V* defined as 

F' ( 9 ) (ipd£) =9(F (V’df )) = J f(x)dx J exp (-j (£, x)) ^ 

= J f(x) exp (-j <£, x)) da; J ip (£) d^ = J f (£) i\) (^) d£ = / (^d^) . 
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We changed the order of integration by means of Fubinis theorem. A similar equa- 
tion holds for F*'. ► 

Example 7. Let q : V — > R be a positive quadratic form: q(x) = 1/2 ^q^XiXj. 
Then 

F(exp(— 2irq(x))dx) = | det q\ exp(— 27rg*(£)) 

where det q = det {qij} and q* = 1/2 QijC& 1S the dual form, i.e., the matrix qj 
is the inverse to q . 

Properties: 

I. The operators F',(F*)' are linear and continuous with respect to the weak 
convergence in 5'(V) and 5"(V'). 

II. They are compatible with the Fourier and Fourier-Plancherel transforms. This 
can be seen from the equation F'[f] = [F(/)] for an arbitrary function / e L\. To 
prove this equation we write 

F '([f])(p) = [/K-F(p)) = J f(x)dx J exp(— j (:r,£)Vd£, 

where we set p £ S. The density on the right side is defined in the product R x R* 
and integrable. Therefore we can change the order of integration: 

F '[f](p) = J f{x)dx J exp(— j <x,£))p(£)d£ 

= J (/ f(x) exp(-j (£, z))da:) p(£) d£ = [F(f)](p). 

III. The operators F', (F*)' are inverse one to another. This follows from Propo- 
sition 1.2. 

The properties of the Fourier-Plancherel transform given in Section 1.4 are fulfilled 
also by the Fourier- Schwartz operators F f and (F*V. We shall write F,F* instead 
of F' and F*' 

Example 8. F'(5o) = [1]. 

Example 9. The function 

x . 1 f 4>(x-y)dy ldy 

iw) = - / — - — = - — (<t> \ x - y )) 

it J y -n y 

is well defined and smooth for any test function <f> and H</> (x) = O at 

infinity. Therefore H (j) belongs to L 2 (M) and the Fourier-Plancherel transform can 
be applied to this function. We can write H (j) = [dx/nx] * (j) where the first factor 
is a distribution. It is easy to check the formula 

F(H(0)(A)dx) = -zsgn(A)0(A), 



(1.5) 




1.4. Homogeneous distributions 
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which means that F'([ dx/x]) = —irisgn£, where sgn£ = ±1 for =t£ > 0. By 
Parseval’s theorem (Section 1.9) this implies that the H, called a Hilbert operator, 
is an isometry in L 2 (M). From (1.5) it follows that H 2 = —I where I is the identity 
operator in L 2 . 

Example 10. We have 

F'((x + (h) _1 da;) = -j£° , F'((x - 0i) _1 dx) = j£° 

where = 1 for ±£ > 0 and = 0 otherwise. It follows from the previous 
calculations and (1.1). 

Example 11. We have for any A G C 

F (H x ) = (j£ + 0)“ A , F (H x ) = ¥W) = (-j£ + 0)“ A - (1.6) 

This implies the group property of Euler kernels with respect to convolution 

H x *H» = H x+ ^. (1.7) 

Example 12. We define the distributions x± k for k = 1, 2, ... as follows 

x± k = x ^ k (*± - (±i) fc_1 r (A + 1) 4 /c_1) ) . 

We have 

F (x; fc dx) = A lim fe (A + k) r (A + 1) [a (£ - o »))- A - 1 - ao fc_1 ' 

a 

= - lim (A + fc) r (A + 1) — (j£ + Of l^fc -1 

A — ► — fc Ofl 

= (fcZiy & + °) fe_1 I ln ( 27r (f - °*)) + OT / 2 ] 
and F ( xZ k dx ) = F (x+ k dx). 

Example 13. The Poisson formula implies F(^25k dx) = By (1.19) for an 

arbitrary h ^ 0, 

F (^2 s khdxj = \h\~ 1 Y^S k/h . 

1.4 Homogeneous distributions 

Let V be a vector space of dimension n. A distribution or a generalized function 
u in V is said to be homogeneous of degree A G C, if it satisfies the equation 
L(e)u = Xu where e = Xid/dxi is the Euler field and L (e) is the corresponding 
Lie derivative (see the next section). This equation means that 

u(L (e) 4>) = —A u(4>) 
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for any test function 0 and any test density p , respectively. This definition agrees 
with the classical one, since for a smooth function or distribution u and test density, 
respectively, function the following equation holds: 

J L(e)u(j) + J uL{e)(j) = J L ( e ) (u<f>) = 0. 

Example 14. The density dx is a homogeneous distribution of degree n and the 
Dirac function 5o((j)dx) = 0(0) is a homogeneous generalized function of degree 
— n, the Dirac distribution <5odx is homogeneous of degree 0. 

Any homogeneous distribution and any homogeneous generalized function is 
tempered. 

Fourier transform 



Proposition 1.4 Let u be a homogeneous distribution in V of degree a ^ 0, —1, 

—2, The Fourier transform u is a homogeneous generalized function on V* of 

degree (3 = —a and can be found by means of the integral 

«(£) = J^T a ((^,x))e H u, (1.8) 

where 

T a (t) = r(a)j _ “(i - Qi)~ a 

and T is an arbitrary cycle in V\ {0} that is homologically equivalent to a sphere 
S mV centred at the origin. 



< Introduce a Euclidean structure in V and take for S the unit sphere. We 



have 



u (4>) = (e H u) |S (ip) , Ip(u) = 




4>(rui)r a 1 dr, 



since u is homogeneous, where |S means restriction to the unit sphere. It follows 
that u £ — » u as e — > 0 in the sense of distributions, where u £ (x) =exp(— e\x\)u(x). 
The distribution u £ decreases fast at infinity and its Fourier transform equals 



Us (0 - (e H u) |s (ipo) , i>o 




exp (— jr (£,w) 



er ) r a 1 dr. 



Set 6 = (£,uj) and use the Euler integral 



f 



exp (— j Or — er) r a 



Mr = T(a)j a (0+je) 



This yields 

U s (0 = r (a) j _a (e d u) (((^, w) + • 
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Take the limit in both sides as e — > 0 and get (1.8) for the cycle T = S. The 
current r a e H u is closed, since we have 



d(r a e H u) = d(e H r a u) = L e (r a u). 



The form r a u is homogeneous of degree 0 hence, the right side vanishes. Therefore 
the integral (1.8) does not depend on the choice of T in the cohomology class of 
the sphere in V\ {0}. ► 

For an arbitrary homogeneous generalized function / in V of degree a > — n, 
the product u = f dx is well defined as a tempered distribution. It is homogeneous 
of degree a + n > 0. The Fourier transform / = F(fdx) is a homogeneous function 
of degree (3 = —a — n and the equation (1.8) reads as follows 

f(0 = T a (£x)f(x)wE, (1-9) 

where c^e = e H dx is the projective volume form. 

Important examples of homogeneous distributions are 

Riesz kernels. Let Q be a positive quadratic form in a space V, say q ( x ) — 
1/2 The quadratic form q* (£) = 1/2 J^q lj ^j (called dual quadratic 

form) is defined in the dual space V*, where the matrix {q tj } is inverse to {%}• 
It is also positive. Consider the family of densities (Riesz kernels) 



R( A) = 



(27 rq) 



X—n/2 



dx 



r(A) 



in V, where dx = dxi A • • • A dx n for a coordinate system x 
This density is locally integrable if Re A > 0. 



(xi,...,x n ) in V. 



Proposition 1.5 The family A i— > R (A) admits analytic continuation on the complex 
plane The distribution R (A) is homogeneous of degree 2 A. 



◄ For a proof we can apply the method of Proposition 1.1 and the identity 



q* (D) q x = A (A — ! + n/2) q X ~\ (1.10) 



where q* (D) = 1/2 ^q lj d 2 /dxidxj. Check the second statement: for Re A > 
n/2 + 1 the density q x ~ n / 2 dx is homogeneous and the equation L(e) (R( A)) = 
2XR (A) holds in the classical sense. This equation holds for all A due to uniqueness 
of analytic continuation. ► 



The function 



Q( A) = 



(2?r q) X ( x ) 
r (A + n/2) 



is locally integrable, if Re A > —n/2 and defines a generalized function in V. 
The family {Q (A)} also admits a holomorphic continuation at C. For any A E C, 
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Q (A) is a homogeneous tempered generalized function of degree 2A, (called a Riesz 
function) . 

Exercise. Check the formula 

S( ~ t)= r(lc + n/2) (~^' (P) ) * = <U ’ 2 '"' 

Define the family distribution R* ( A) = (2 / 7rq*) X ~ n/ ' 2 dx/T (A) and of generalized 
functions Q (A) = (2nq*) X /V (A + nj 2) in the dual space V* by means of the dual 
quadratic form q* . Take analytic continuation of these families at C. 

Proposition 1.6 For any A G C the Fourier- Schwartz transform of the Riesz kernel 
of degree A is the dual Riesz function of degree —A and vice versa: 

F {R (A)) = Q* (—A) , F (Q (A)) - R* (-A) . (1.11) 

In particular, 

f (2Trq) x ~ n/2 dx \ = (2nq* ($))~ A 

V r W ) r(n/2-A) 

for all A such that nj 2 > Re A > 0. 

◄ Introduce the Euclidean structure in V such that q (x) = \x\ 2 ; then we have 
q* (£) = |£| 2 . Apply Proposition 1.4 to the homogeneous distribution u of degree 
a = 2 A such that u\S = 1. Then we only need to evaluate the integral 

h(0= [ ((t;,u) + oi) 2X dLu, 

J S 

where do; = e H dx is the Euclidean area on the unit sphere S. The integral is 
spherically symmetric and reduces to a Beta function: 

I (0 = 5 n _! (1 + exp (jA)) 2 ~ X B (1/2 - A ,(n- 1) /2) |£f 2A > 

where s n _i is the area of the unit sphere in E n_1 . This implies the first equation 
(1.11). The second equation follows from the first one. ► 

Powers of arbitrary quadratic form. Let q be an arbitrary real non- 
singular quadratic form in V. One can define the family of locally integrable 
functions Q+ = q\/T (A + 1) T (A + n/2) where A G C, Re A > —1. This fam- 
ily can be extended to the whole plane C by means of the identity (1.10) which 
holds for any quadratic form q. This implies that the family of functions q+ has 
meromorphic continuation with occasional poles at the points A = — 1, — 2, . . . and 

A = —n/2,— n/2 — 1, Another method is similar to that of Example 5. We 

define powers of q as follows: 

(q ± 0z) A (p) = Urn f (q± ei) X p, p € K (V) . 

E^+0J 



( 1 . 12 ) 
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The power ( q ± ei) X is defined by means of the branch of the power function 
£ A = exp (A In C), where ln£ is the branch of the logarithm in C\ (— R+) that is 
real for positive £. 

Proposition 1.7 The limits (1.12) exist for any p and define holomorphic families 
A i-> (q ± (h) A A G C, of generalized functions in V. 

◄ Choose a positive quadratic form p = p(x) in V Fix e > 0 and consider the 
family of functions g A where q £ (x) = q (x)+ep. It is holomorphic for A £ C since q £ 
does not vanish in V. It fulfils the identity q * ( D ) q x (x) = A (A — 1 + n/2) g A_1 (x) 
where q* is the homogeneous second-order operator dual to the form q £ . Therefore 
we have identity 

The right-hand side has meromorphic continuation to the half-plane {A;ReA> 
— 1} with occasional pole A = 0, since the dominator vanishes. Integrating by 
parts yields that the numerator also vanishes as A = 0. Therefore there is no 
pole at this point and we get analytic extension of both sides to the half-plane 
{Re A > —1}. Using again (1.13) gives analytic continuation at {A; Re A > —2} and 
so on. Another occasional pole isA = l — n/2. Integrating by parts yields again 

l ql~ n > 2 (q* (D) p) = 0, 

which means that there is no pole at all. Now we take £ — > 0; the dual form q * 
tends to the dual form q * since q is non-singular. The right-hand side of (1.13) 
has a limit for Re A > —1 hence, the left side also does. Repeating the above 
arguments, shows that the family q * has a limit as e — ► 0. It coincides with the 
function (q + Ck) A for Re A > 0. This implies that the family (q 4- 0z) A has analytic 
continuation on the whole plane. The same is true for (q — Ck) A . ► 

Comparing the power functions, we get (q =b 0z) A = g^+exp (±7rAz) g* , where 

= (-«)+• 

For an arbitrary natural integer k we define the generalized function 
[<T k ]= \ {q + 0t)~ k + {q - 0i)~ k . 

Note the equation 

= { ^j^-[(q-0i)- k -(q-0i)- k ' . 



(q) 



(1.14) 
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1.5 Manifolds and differential forms 



Differential forms and orientation. Let V be a vector space and x = 
(xi,...,x n ) be a coordinate system in V. The differentials daq, . . . , dx n are dif- 
ferential forms of degree 1 in V. For an open set U C V and consider the algebra 
O ( U ) of infinitely differentiable functions in U. Take the exterior algebra Q ( U ) 
generated by the forms daq, . . . ,dx n over O ([/); it is called the algebra of even 
differential forms in U. Any even form a of degree k is an element of 0 ({/); it can 
be written as 

ot — (x)dx il A • • • A dx ik , (1.15) 



where A is the symbol of the exterior product. This operation is bilinear and 
fulfils the condition (3 A a = ±a A (3 where the sign is — if both forms a , (3 have 
odd degrees and + otherwise. The coefficients belong to 0 ({/); they are 

uniquely defined if the sum (1.15) is taken only over the set of indices such that 
i\ < • • • < ik^ To write an even form in another coordinate system y -- (yi, ... , y n ), 
we need only to replace each generator of the algebra Q(U) by means of the 
standard rule 



dx, 



dXi = X df dyj ' * = !> • • ■ * 



(1.16) 



Definition. Let X be a smooth manifold; the structure of the manifold is given by 
a set (atlas) of maps (U,x), where U is an open set in X and x =- (aq, . . . ,x n ) 
is a coordinate system in U. An even differential form in X is a form oljj defined 
for each map (U,x) in X such that for any two maps (U,x) , (V,y) the forms au 
and ay are related in U D V by the rule (1.16). Let Q ( X ) be the space of even 
differential forms in X. The exterior differential of the form a of degree k is the 
globally defined form da of degree k + 1, defined locally as 



da = da^ 



.,i fc A dx i± A • • • A dx ik . 



Definition. An orientation on a manifold X is a choice of an atlas of maps such 
that the Jacobian det dy/dx is always positive for any two maps (U,,x ) , (V,y) in 
the atlas. If X is an oriented manifold of dimension n, an orientation form on X 
is an arbitrary even n-form a on X such that a (d/dxi , . . . , d/dx n ) > 0 for any 
coordinate system aq, . . . , x n that belongs to the atlas of the orientation. 

Example 15. Let V be a vector space of dimension n. The set of one-dimensional 
subspaces L C V has the structure of an (algebraic) manifold of dimension n — 1 ; 
it will be denoted P (V) (the projective space ofV). Take a system of coordinates 
xi , . . . ,x n in V; the set Uj of lines L such that Xj ^ 0 in L is an open set in 
P (V). The coordinates x \, . . . , Xj_i, #j+i, . . . , x n of the intersection point y ( L ) = 
L n {xj = 1} form a local coordinate system in Uj,j = 1, . . . ,n. The manifold 
P (V) is orientable if n is even and is not orientable if n is odd. 
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Integral. Let A be an oriented manifold; the integral f x a is well defined for 
any even differential form on X with compact support and continuous coefficients. 
For any even form a the “Stokes” integral formula holds, 




where dU is the smooth boundary of a compact oriented open set U C X. The 
orientation of dU is defined by means of the form df where / is a function that 
fulfils the conditions: / < 0 in U and / = 0, df ^ 0 on dU. This means that a local 
coordinate system y = (y 2 , • • • , y n ) belongs to the atlas of the orientation of U, if 
the coordinate system y\ = / (x) , Z 2 , . . . , z n belongs to the atlas of orientation of 
U and yj = Zj\dU, j = 2, . . . , n. 

Definition. The form a is called odd , if there is one more factor sgndet ch//<9x in 
the transformation formula from coordinates y to coordinates x. Any odd form 
with compact support can be integrated over the manifold X with no orientation. 
A volume form a in X is an odd form such that a = adxi A • • • A dx n and a ^ 0 
for any local coordinate system. There is no Stokes theorem for odd forms. See 
[18] for more information. 

Example 16. A measure, a charge, a density, a surface density = area element, a 
curve element are all odd forms. 

Critical points. Let / : X —> Y be a smooth mapping of manifolds of local 
dimensions dim (A, x ) , dim (Y, y). A point Xo E X is called critical , if the rank of 
the linear mapping df (xo) is l ess than dim (Y, f (xo)) Note that the mapping df 
(xq) is given by the Jacobian matrix {dyj/dxi} in local coordinate system {xi} 
at xo and {yj} at / (xo) respectively, and rankd/ (xo) = rank {dyj/dxi}. A point 
yo G Y is called a critical value of / if yo = / (xo) for, at legist one critical point 
xo- By Sard’s theorem the set of critical values has zero measure in Y. 

Degree of a mapping. A mapping / : X —> Y of manifolds (or topological 
spaces) is called proper , if for any compact set K in Y the set / -1 ( K ) is compact in 
X. Let / be a proper mapping of smooth oriented manifolds of the same dimension 
n = dim A = dimY. Take an arbitrary non-critical point z G Y and consider the 
sum 

deg(/,z) = 53sgn^| /(x)=z . (1.17) 

The sum is finite since the mapping / is proper and each sign is well defined since 
z has non-critical value. 

Proposition 1.8 If Y is connected, the number deg ( f,z ) is constant. 

M If we move the point z along a generic curve in Y, then either each term 
in (1.17) stays constant or two terms 1 and —1 disappear or appear simultan- 
eously. ► 

The number (1.17) is called the degree of the mapping / and will be denoted 
by deg (/). The degree depends on the orientations and changes its sign if an 
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orientation is changed to the opposite one. The basic property of the degree is 
as follows: if a is an even differential n-form on Y with compact support (more 
generally, an integrable form on 7), then 

/ /* (<*) = deg (/) f a. 

Jx Jy 

Example 17. Consider the mapping p : C — > C defined by a complex polynomial 
P — P( z )- Critical points of p are the roots of p'. The degree of this mapping 
coincides with the degree of the polynomial p. 

Tangent vectors and covectors. Let x be a point in a smooth manifold 
X , O x be the algebra of smooth functions defined in a neighborhood of x. Any 
functional t : O x — > R satisfying the equation t ( ab ) = t(a)b ( x ) + a (x) t ( b ) , a, b G 
0 X is called a tangent vector in X at x. The space of all tangent vectors at x is 
denoted T x (tangent space at x ) . The dual vector space T* is called the cotangent 
space at x. An element r of the cotangent space is called a differential form of 
degree 1 or covector at x. Let Y be a submanifold of X. A vector t is tangent to Y 
if t (a) = 0 for any function a G O x that vanishes in Y. Any functional t G T * such 
that t (t) — 0 for any tangent vector t to Y is called a conormal covector to Y at 
x. The union of spaces T (X) = U T x , T* (X) = U T* are called the tangent bundle 
and the cotangent bundle, respectively. If X is an open set in a space V = M n , 
the tangent bundle is isomorphic to X x V and the cotangent bundle has the 
structure T* (X) = X x V* where V* is the space dual to V. Let O (X) be the 
algebra of smooth functions defined in X. A linear operator t in this algebra is 
called a tangent field , if it satisfies the Leibniz condition: t (ab) = t(a)b + at (b). 
Any tangent field t defines at each point x G X the tangent vector t x : O x — > R 
such that t x (a) = t(a) for functions a defined on X. A differential form a of 
degree A; at a point x G X is a mapping T x x • • • x T x — > C, i.e., a (fy, . . . , t n ) G C 
that is linear in each argument and changes its sign if two neighboring arguments 
permute. If the form a is given by (1.15), then 

a(ti,...,t k ) = a h (si) A • •• Adx ik (s k ) , (1.18) 

where the inner sum is taken over all permutations of the vectors 

ti, . . . , tk and 7 r is the parity of the permutation. Vice versa, any differential form 
at a point x can be written as the sum (1.15) with some coefficients a... G C. A 
differential form a of degree k on the manifold X is an arbitrary smooth family of 
differential forms a x , x G X of degree fc; the condition of smoothness means that 
the coefficients a... in (1.18) are smooth functions of x in any chart. 
Contraction. Let s be a tangent field in X and a be a differential form of degree 
k > 0. The contraction of a by means of s is the differential form (3 or degree k — 1, 
denoted by s b a, such that 

P if 1? • • • j tk~ l) — Qf(s,fy,..., tk— l) • 

We have always s h {s b a) = 0 and t\ b (£2 L • . . (tk l - a)) = a (fy, . . . , tk ). 
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Example 18 . Let xi,... ,x n be a coordinate system in a space V and e = J2 %id/dxi 
be the Euler tangent field. Then 

e H dxi A • • • A dx n = (— l ) z_1 Xidxi A • • • A dxi-i A dxi+\ A • • • A dx n . 

Given a tangent field v, there is defined the Lie derivative L (e) along v of an 
arbitrary even differential form a, 

L (v) a = d (v h a) + v h da. 

The operator L (v) preserves the degree of the form and commutes with the exterior 
differential. For a function a (form of zero degree) the Lie derivative coincides with 
the action of the tangent field: L (v) a = v (a). 

Quotient of forms. A form (3 of degree k is called non-degenerated at x if 
(3 (£i, . . . , tk) ^ 0 for some tangent vectors at x. Let (3 be a non-degenerated form 
and a be an arbitrary differential form of degree n = dim A. The Leray form or 
quotient a/ (3 is an arbitrary n — k - form 7 such that (3 A 7 = a. Let / 1 , . . . , /*. 
be some smooth functions in X such that the form /? = d/i A • • • A dfk is non- 
degenerated. Then the quotient a/f3 is uniquely defined up an additive term of 
the form d/iAftH b dfk A /?&. 

Definition. A manifold X is called a Riemannian manifold if there is given a 
smooth positive symmetric bilinear form g (called a metric) on the tangent bundle 
T (X). These conditions mean that for any tangent vectors u,v at a point x we 
have g(v,u ) = g(u,v), g(u,u) > 0 , if u ^ 0 , and for any coordinate system 
x = (xi , . . . , x n ) in an open set U C X the coefficients gij (x) in the local expression 
g = J2 9ij ( x ) dxidxj are smooth functions in U. For any C 1 -curve 7 given by the 
parametric equation x = x (t) , 0 < t < 1, the length I 7 I is defined to be 

hl = l 9' /2 («><i‘=E^ 

For any points y,z £ X the distance dist (y, z) is the length of the shortest curve 
that joins these points. A geodesic curve in A is a curve 7 which is locally a 
shortest way between any two points 2 /, z E 7 whose distance is sufficiently small. 

The Riemannian metric g gives sense to the length of a vector, to the angle 
between two vectors and induces the area density on any smooth submanifold 
Y C X. A submanifold Y is called totally geodesic, if any geodesic curve that is 
tangent to Y is contained in Y. 

See more information in [98] or other textbooks on differential geometry. 

1.6 Push down and pull back 

Distributions and generalized functions have different behavior under change of 
variables. Take a smooth invertible mapping / : U — > V defined in an open set 
U C W. Then any test function ^ E D (V) is transformed to a test function 
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4> £ D ( U ) by 0(x) = /* (^) (x) = if (/ (x)). The pull back (inverse image) of a 
test density p — if dx is equal to cr = f* ( p ) = (f |det dyjdx | dx where dy/dx is the 
Jacobian of the mapping /. The push down (direct image) of a distribution u and 
pull back (inverse image) of a generalized function v are defined in the natural 
way: for u G D' ( U ) , v E K f ( U ) we set 

/* W) WO = u (/* WO) , /* W) W) = v (9 * (p)) » where g = / _1 . (1.19) 

If v = v (x) is a locally integrable function, we have /* ( v ) (x) — v (f (x)). 

Example 19. Take a linear invertible transformation A : M n — » M n . For the 
Dirac distribution Jodx we have A* (Sodx) = Sody, unlike the equation A * (So) = 
|det A] -1 Jo for the Dirac function. 

These constructions are generalized for more general mappings of smooth 
manifolds and distributions and generalized functions defined on such manifolds. In 
the general case, a generalized function may have pull back, whereas a distribution 
may have pull down. We consider here only the first operation. 

Pull down of a test density. Let X, Y be smooth oriented manifolds and 
/ : X — > Y be a submersion, i.e., a smooth mapping without critical points, i.e., 
such that rankd/(x) = dim (F, /(x)) for any x G X. Take an orientation form 
a in X, an orientation form p in Y and consider the pull back /* (p); this is a 
non-degenerated form in X of degree dimF. Take a test density ip G K (X) and 
consider the quotient <p/f* (p); it is a differential form of degree dimX — dimX. 
Each fibre X y = / -1 (y) is a smooth manifold of dimension dimX — dimX, since 
/ is a submersion. Define the orientation of X y by means of the volume form 
<r/f* (p) and consider the integral 



^ = J x .f%ry 

The integral converges, since the support of (p is a compact set. The function if is 
smooth in Y, since / is submersion. The product /* (<p) = if (y) p is defined in Y 
and does not depend on the choice of p (but depends on orientations): for another 
choice p' we have p' = ap where a > 0 and (p/p' = a _1 (p/p, if' = a~ x if, if'p' = ifp. 
The form /* (<p) is called the push down of <p. We can say in a non-formal way: push 
down of a form is the result of integration of the form along fibres of the mapping. 
It is easy to check that the form /* (<p) is smooth and has compact support, i.e., 
f m (< p )€ K ( Y ). 

Pull back of a generalized function. For an arbitrary generalized function 
v in Y we set 

/* W) (</>) = v (/* (< p )) . 

The functional u = /* (v) is well defined and continuous on the space K (X), i.e., 
u is a generalized function in X. It is called the pull back of v. 
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Example 20. Let V be a vector space and U be an open subset of V. Take a 
submersion / : X — ► M and the Dirac function v = Jo- By the above definition the 
pull back is denoted 5 (/) = /* (So), is well defined and 



<5 (/)(</>) 



f & = lim 1 [ v . 

Jf = o d/ £->+o 2e J\f\< £ 



Here t = f (x) is the coordinate form of the mapping /, dt is the orientation form 
in R. Formally d / = /* (dt). If X is an open set in a Euclidean space E, dV is the 
volume form in E and p = adV, then we can write 

‘w-Lm* 

where d5 is the hypersurface area form in E. If we replace / by the function 
g = af , where a ^ 0 is a smooth function, then 

S(g) = ^-S(f). 

\a\ 



Pull backs of derivatives of the Dirac function can be written in a similar way: 
taking v = 5^ gives the equation 



s {k) (/) = 




\t = 0 



which agrees with (1.14). Replacing / by g = af for a constant a / 0, gives 

S (k) (9) = jXk^ (/). 

\a\ a K 



Example 21. Let / be a real smooth function in an open set U C V such that 
df 0 as / = 0. We define the generalized Euler kernel generated by / as 

H f {p) = T(X}~ {P) = f(A) J f> o /A_ ^ dx ’ D(U), Re X> 0. 

In the case n = l,f = x this is just the Euler kernel as in Section 1.2. It possesses 
similar properties: 



Proposition 1.9 The family Hj has analytic continuation on C with values in 
K'{U). The functional gW (/) = HJ k is supported by the hypersurface Z = 
{/ = 0} for k = 0,1,2, 

◄ Define the function 
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where the hypersurface is oriented by the form df. The function r has compact 
support, since p has such support. It belongs to C°° in a neighborhood of the 
point t = 0, since df ^ 0. Therefore we can apply the Euler distribution to the 
test function r. By Fubini’s theorem, 

H x (r) = ReA > 0. 



The left side has analytic continuation in the complex plane. This gives analytic 
continuation for Hj. We have 



Hj k (<j>) = <5 (fe) (r) = r {k) (0) 




(j)dx 

~df 



► 



1.7 More on the Fourier transform 

In this section, knowledge of rudiments of Lebesgue integration theory is assumed. 
Let X be a space supplied with a Lebesgue measure dx . A (measurable) function 
/ : X — » C is called Lebesgue integrable in (X, dx), if the Lebesgue integral 

ll/lli = / l/|dx 

is finite. A function / is called negligible if / = 0 almost everywhere; this condition 
is equivalent to the equation \\f\\i = 0. Any equation like / = g in Lebesgue’s 
theory means only that f — g is a negligible function. 

For an arbitrary number p > 1 the notation L p = L p (X),p > 1 is used 
for the C- vector space of functions in X such that the function \f\ p is Lebesgue 
integrable. More accurately, this is the quotient space of such functions modulo 
the subspace of negligible functions. This is a Banach space with the norm 

ii/ii P = (/ i/rdz) 1/P . 

In particular, the space of square- integrable functions in X, is a Hilbert 

space with the scalar (inner) product 

(f,g) = J fgdx. 

Remind the inequality: |(/, g)\ 2 < ||/|| 2 ||g|| 2 , which holds for any scalar product. 
It is equivalent to the triangle inequality \\f + g\[ 2 < (ll/ll + IMI) 2 for L 2 -norms. 
We need few more facts from the Lebesgue theory: 

Theorem 1.10 (Dominated convergence theorem) Let F be an integrable function 
in V and fi,i = 1,2, . . . a sequence of (measurable) functions such that \fi\ < F 
and fi—>f almost everywhere in X. Then 

[ fidx -> [ f dx. 

Jv Jv 
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Theorem 1.11 (Fubini theorem) Let X , Y be spaces endowed with the Lebesgue 
measures dx,dy, respectively, and f be a function in X x Y integrable with re- 
spect to the measure dxdy in X xY. Then the function f(’,y) is integrable in X 
for almost all y , the function g (y) = f x f (x,y) dx is integrable in Y and 

dy f(x,y)dx = f (x, y) dxdy. 

JY JX JXxY 

From Fubini’s theorem we conclude that 

/ ^ j f(x,y)dxjdy = J ( j f(x,y)dy S j dz, 

i.e., we may change the order of integration for any integrable / in X x Y. 

Let V be a finite-dimensional vector space. Fix a system of coordinates and 
consider the corresponding volume form dx in V. Lebesgue theory is then applica- 
ble in V. We say that a function g in V is a step function, if it is equal to a linear 
combination of indicators of cubes Q C V. Recall that the indicator function x of 
a set G is equal to 1 in G and 0 otherwise. 

Theorem 1.12 (Density theorem) For any p> 1 the set of step functions is dense 
in the space L p (V). 

The Fourier transform of a function / G L\ (V) is defined as in Section 1.1. 

Theorem 1.13 (Parseval) For an arbitrary function /gLi (V)flZ/ 2 (V) the Fourier 
image f belongs to L 2 (V*). Then the equations 



ll/ll 2 = ll/ll 2 , 


(1.20) 


(f, 9 ) = (/> 9 ) 


(1.21) 



hold for arbitrary g G L\ D 1/2- 

In other words, the operator F : / i-> / is an isometry in L 2 (V) —> L 2 (V*). 

◄ Consider the case V = M for simplicity. Take a finite interval, say I = 
[a, b] C R and denote by hi the indicator of this interval. First we check (1.21) for 
f = hi,g = hj for arbitrary intervals 7, J. We have 



hi 



-f 



exp(— j£x)dx 



jf 



[exp (— j&£) — exp (-ja£)] 




a 4" b 



sin(7r£(6 — a)) 
7r£ 



2 
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We see that the right side is a holomorphic function in the whole plane C (more- 
over, a bandlimited function, see Section 1.8). Similarly for J = [c, d] we have 






J hj hjd£ = J 



-j^ _ e -jc£j ? 

e j (d-b)€ _ e j (d-a)€ _ e ) (c-6)£ _|_ g j (c-a)£ 



(2tt0 2 ‘ 



There is no pole at the point £ = 0 and we can integrate over the line £ = rj — i. 
The first term in the bracket gives zero after integration, if d — b < 0, since the 
function exp (j (d — b ) £) decreases rapidly in the bottom half-plane. Otherwise the 
integral of this term is equal to j 2 (d — b) (apply the residue theory). Therefore the 
result of integration equals 



J hi hjd £ = - (d - 6 ) + + (d — a) + + (c — 6 ) + — (c — a) + , 

where we set e + = e if e > 0 and e+ = 0 otherwise. It is easy to see that the right 
side coincides with the length of the interval /fl J which is equal to f hjhjdx. 
This proves (1.21) for f = hi,g = hj and for arbitrary step functions / and g. 

Take an arbitrary / G L\ fl L<i and an arbitrary t > 0; set /*(#) = f(x) 
for \x\ < t and f t (x ) = 0 otherwise. We have ||/ — ftlh — > 0 and ||/ — f t || 2 — > 0 
as t oo. By Theorem 1.12 we can choose a function h t that is equal to a 
linear combination of indicator functions of some intervals [a, 6] C [—t,t] such 
that || ft — h t || 2 < l/t. We have ||/ t — htlh — (2t) 1//2 /t = ( 2/t ) 1//2 by the triangle 
inequality, consequently ||/ — h t \\ 1 — > 0 as t — > oo. By 1.4 ^ —> f uniformly. On the 
other hand, \\f — h t || 2 — > 0, which yields that {h t } is a Cauchy sequence in L 2 (K). 

By (1.20) is a Cauchy sequence in the space L 2 (M*) and by completeness of 

this space h t (f) in L 2 (M*). We conclude that </> = /, since the same sequence 
converges to / uniformly. Finally 



; 



2 = W 2 = lim 



= lim \\h t 



2 j 



which proves (1.20). The equation (1.21) follows from the identity 2Re(A/,g) = 
|| A/ + g|| 2 — |A| 2 1|/|| 2 — ||g|| 2 , where A is an arbitrary complex number. ► 



Theorem 1.14 (Plancherel) For an arbitrary function f G L 2 (V) the sequence of 
functions 

Fr(f)(0 = [ exp(— j£x)/(x)dx, r > 0 

JB(t) 

converges in L 2 (V) to a function f as r ^ 00 where B (r) means the ball of radius 
r. The limit satisfies (1.20). 
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◄ Let again n — 1. The integral F r (/) equals the Fourier transform of the 

truncation f r of /. The sequence {f r } converges in L2 to / as r — > 00, hence it 
is a Cauchy sequence in L2. We can apply the Parseval equality, since f r E L\\ 
II /r — /sib = || fr — /sib- It follows that / r ,r = 1, 2 , ... is a Cauchy sequence L2. 
Denote by / its limit. We have || f r — f\\ — > 0, which yields ||/|| 2 = ||/|| 2 . ► 

We call the function / the Fourier- Plancherel transform of /. This transform 
is a continuation of the Fourier transformation, since the function / coincides with 
(1.4) if / E Z/inZ/2. The continuation is unique and keeps the Parseval equality. We 
use henceforth the same notation F : f f for the Fourier-Plancherel transform. 
The Parseval equality means that this operator is unitary as an operator in L2. 
The adjoint Fourier transform F* possesses similar properties and can be extended 
to the space L2 (R*). 

Proposition 1.15 For any ip E Z/2(M), ^ £ L 2 OR*) we have the equation 

= ( 1 . 22 ) 

◄ First assume that both functions ip, are integrable and change the order of 
integration as follows 



= J J exp(-j (t,x))<p(x)dxip(t)d£ 

= J J ¥’(a;)exp(j (^,x))V>(0<^da: = {<p, 

For arbitrary square-integrable functions we can apply this equation to 
the truncated functions ip r ^r • Then we pass on to the limit as r — ► 00 and 
get (1.22). ► 

Theorem 1.16 The operators F, F* are mutually inverse, i.e., 

F*F = I, FF* = I, 

where I means the identity operator in Z/ 2 (V) omd in L 2 (V*). 

◄ The second equation is similar to the first one. The composition F*F is an 
isometry. Therefore it is sufficient to prove the first equation on a dense subset 
of Z/2. For indicator functions we can change the order of integration and apply 
Parseval’s theorem: 



(F*F (hj) , hj) = (F (hj) , F (hj)) = ( hj , hj ) . 



This implies that F*F(h/) = hi and by continuity of the operator F*F this 
equation holds for any function / E L2 (R*). ► 
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Convolution. The integral 

(/ * a)(x) = Jf(x - y)g(y)dy (i-23) 

is called convolution of integrable functions /, g. The mapping (/, g) i— > / * g is a 
bilinear, commutative and associative operation. It satisfies 

\f*g\<\\fh\\gh (1-24) 

which follows from the triangle inequality. 

Proposition 1.17 If f,g E L\ the integral (1.23) converges almost everywhere and 

II/*5||i<II/IIiNIi. 



◄ By Fubini’s theorem and by changing the variables z — x — y we get 

III f(y)g(x - y)dy^dx < j \f{y)\\g{x - y)\dxdy = J \f(y)\dy J \g(z)\dz 



Theorem 1.18 For arbitrary f,geLi we have 

F(f * g) = F(f)F(g). (1.25) 

If moreover, /, g E L 2 the symmetric equation is also valid: 

F(fg) = F(f)*F(g). (1.26) 

These equations hold for the conjugated transform F* as well. 

M First we prove (1.25): 

F(f * g)(0 = J exp(-j (£,x)) J f(x - y)g(y)dydx 

= J exp(— j (£,*)) J exp(— j (£,y))f(z)g(y)dydz 

= J exp(-j (£,z})f(z)dz J exp(— j {£,y))g(y)dy. 

The coordinate change z = x — y is eligible, since the integrand belongs to 
the Lebesgue space Zu(V x V). The right side is equal to F(f)F(g). To prove 
(1.26) we apply (1.25) to the adjoint Fourier transforms of the truncated func- 
tions F(f) r , F(g) r . These functions are integrable since they are bounded; this 
yields 



F*(F(f) r * F(g) r ) = F*(F(f) r ) F*(F(g) r ). 



(1.27) 
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By Plancherel’s theorem the sequence F*(F(f) r ) converges in mean to F*(F(f)) 
as r — > oo and F*(F(f)) = f by Theorem 1.16. Similarly F*(F(g) r ) — > g in mean. 
The right side of (1.27) converges to the function fg in the space L\ by the triangle 
inequality. Apply the Fourier transform to both sides and get uniform convergence 
of the sequence FF*(F(f) r * F(g) r ) — * f g. On the other hand the convolution 
F(f) r * F(g) r has compact support and is bounded by virtue of (1.24). It belongs 
to L 2 hence FF*(F(f) r * F{g) r ) = F(f) r * F(g) r by Theorem 1.16. By (1.24) 
F(f) r * F(g) r converges uniformly to F(f) * F(g) which implies (1.26). ► 

Convolution of distributions. The convolution can be correctly defined for 
some pairs of distributions u,v in V. Let </> be a test function; one set u * </> (x) = 
u where (*) = 4> (x — ■). The function u*(f> is always smooth but need not 
have compact support. In general the support of the convolution u*(f) is contained 
in the sum of supports of <f> and u. The formula 

U*v(4>) — u (v * 0) 

makes sense if v has compact support. The convolution u * v is well defined as a 
distribution in V. If u has compact support, the convolution u*v can be defined by 
a similar formula. If both u , v have compact support, the two definitions coincide. 
Another option is: V = M and the supports of u and v are contained in a ray 
4- a for some real a. Then the support of u * v is contained in M+ + 2 a. 

Example 22. We have supp H x C M+ for the Euler kernels, hence the convolution 
of kernels H x and H M is always defined. 

Proposition 1.19 H x * H ^ — JT A+M for A, g G C. 

◄ For Re A > 0, Re g > 0 this is just the Euler identity for Gamma and 
Beta functions. It remains true for all A, //, due to uniqueness of analytic con- 
tinuation. ► 

Example 23. Equation (1.25) holds also for some pairs of distributions f,g such 
that the convolution / * g is well defined. In particular, formula (1.5) implies that 
(1.25) holds for the distribution / = [dx/x] and any density g E T 2 (R). 

Poisson formula 

Theorem 1.20 Let f be a function in R fulfilling the inequalities 

\f(x)\<C(\x\ + l)~ a and |/(0|<C'(|£| + ir° 
for some a > 1,C > 0. Then the equation holds, 

z z 



(1.28) 
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◄ Note that the inequality for / implies absolute convergence of the right side 
of (1.28) and integrability of / on the whole line. Therefore the Fourier transform 
of / is well defined. The inequality for / implies absolute convergence of the left 
side and continuity of /. The series 

9(x) = ^2f( x + k) (1.29) 

kez 

also absolutely converges and the sum g is a continuous function in [0, 1 ]. Take the 
Fourier series 

9 ( x ) = ^2 9n exp (j nx) ,g n = g(x) exp (-j nx) d x = f (x) exp (-jnx) dx. 

nd-lf " 0 v M 



Substitute the series (1.29) in the integral and find that g n = f (n) .. This implies 
that the series \g n \ converges and the first equation (1.29) holds at each point. 
Taking x = 0 yields (1.28). ► 

Remark. The conditions /, x 2 /, /" G L 2 (M) also imply (1.28). 



1.8 Bandlimited functions and interpolation 

Theorem 1.21 (Paley- Wiener) If the support of a function f G Z/2(M) * s contained 
in the interval [a, b] for some a, b, then its Fourier transform p = F (/) has analytic 
continuation in the complex plane that satisfies the inequality 

|y>(C)| < Cexp(27rmax(6ry, arj )) , ( = £ + irj. (1.30) 

If a function p G L 2 (R) has a holomorphic continuation in C that fulfils (1.30), 
then the function f = F* (p) vanishes almost everywhere in M\ [a, b], 

◄ The kernel of the Fourier integral 

f(0= [ ex P / (x) &X 

J —a 

has holomorphic continuation exp (— jOr) which is bounded in the segment [—a, a]. 
We have |exp(j£x)| = exp (fZnrjx) < exp (27rmax ( 677 , arj)) which proves the first 
statement. 

Suppose that p fulfils (1.30) and f = F* (p). We show now that (f,g) = 0 
for an arbitrary function g G D (R) with compact support supp g C [c, d] where 
c > b or d < a. This will imply the second statement. By Parseval’s theorem, 
(/, g) = (p, V’) where ^ — F {g). By the first statement the function has analytic 
continuation in the complex plane which fulfils an estimate like (1.30). Suppose 
that c > b and consider the function V 7 * (C) — V 7 (C); it is analytic too and fulfils 

1^(01 < C'OCI + l) _2 exp(-27TC?7), r\ > 0. 
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By Cauchy’s theorem we can move the chain of integration: 

{<p, 4>) = J v (0 (0 = J v (f + vn) i’* (£ + l v) <%■ 

According to (1.30) we have for an arbitrary rj > 0, 

1^(0^* (01 < C (|C| 2 + l) exp (27r (by — cr })) , 

which implies \((f,ip)\ < C exp (2n (b — c) rj). Taking 77 — > 00 , yields {<p,ip) = 0, 
q.e.d. In the case d < a we take rj — > — 00 . ► 

Example 24. Consider the function ip = (1 — x 2 )^ -1 / 2 /r (A + 1/2), where A E C. 
It belongs to Li (R) for Re A > —1/2. Take the Fourier transform 

i> = n X J ^ p ^ X + exp (-j£r) da: = |^|“ A Ja(0- 

The factor J\ is called the Bessel function of order A. By Example 4, the left-hand 
side has analytic continuation at C, hence the Bessel function A t— ► J\ has analytic 
continuation too. 

Definition. A function / E Z/ 2 (R) is called an a-bandlimited function for some 
number a > 0, if supp / C [—a, a]. 

Example 25. The function sincx = sinx/x is a-bandlimited with a = 1/27 r. 

Properties: 

1. If g e I /2 and supp <7 C [ 6 , c] for some b < c, then the function /(x) = 
exp (— (b + c)7rix)g(x) is a-bandlimited for a = (c — 6)/2. 

2. If / is an a-bandlimited function and g is a 6 -bandlimited function, then fg 
is an a + 6 -bandlimited function. 

3. If / is an a-bandlimited function, then shift Thf is again an a-bandlimited 
function. The derivative /' is a-bandlimited too and \\f '\\2 < 27ra || /* || 2 - For 
a proof we take in account the equation f r = j £/. From this \f'\ < 27ra|/|, 
since / vanishes for |£| > a and by Parseval’s theorem 

Wf'Wl = J |/fd^<(27ra) 2 | |/| 2 dC = (27ra) 2 ||/|| 2 . 

4. Any bandlimited function can be interpolated in an arbitrary finite interval 
from its complement: 



Proposition 1.22 For arbitrary positive a and b any a-bandlimited function (p E 
L 2 (R) can be interpolated in (—6, b) as 



<£(C) = exp ( 27 ray/b 2 - C 2 ) J - J 



-b sin ^2na^y 2 — 6 2 ) 



n(y- 0 



0(y)dy (1.31) 



where the branch is defined by Re \Jh 2 — Cfi > 0. 
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◄ Denote T + = [h, oo),r_ = (— oo,6],r = T + U (— T_). Take the meromorphic 
form 

—2 irayjb 2 — z 2 ) (z) dz 

i( z ~0 

defined in C\T. It tends to zero at infinity due to the Paley- Wiener theorem and 
is absolutely integrable on each side of T. Take the e - neighborhood T± (e) of r± 
and the chain 7 (e) = 9r_|_(e) — <9r_ (e). By the residue theorem for a point 

C e c\r + (£)ur_ ( e ), 

exp (- 7 ry/b 2 - C 2 ) (j) (C) = - f a. 

Choose now ( G T. For two close points (± G 7 (e) such that d=Im£± > 0, 
Re >0, we have 




exp (-TT^-Ci) — exp )*2t sin 6 2 ) , 

where 2£ = £+ + C- • 

In the case Re£± <0 we get the quantity « — 2zsin ^7r^/ £ 2 — b 2 ^j instead. 



Therefore 



as c — * 0. 



7 *7 

7 7 ( e ) Jr 



sin ( 7 r\/z 2 — 6 2 ) 

tt(2:-C) 



0(z)dz 



Corollary 1.23 If (p is an a-bandlimited function, then for any b > 0 we have 

f exp (—27 xa\/b 2 — x 2> ) <p(x) dx < f \<j) (x) | 2 dx. 

7-6 V 7 7|x|>6 

◄ Write the equation (1.31) as 



/ , X f sin (tt^/j/ 2 - bA 

exp(— 27 t ay b 2 — x 2 ) (j) (x) = / — — - 

\ J J r 7r(y-x) 

f i>( y)dy 
J Tr(x-y)’ 



<t> (y) d y 



\x\ < 6, 



where xp (y) = =Fsin (ny/y 2 — b 2 ^j <j> (y) for y G T± and xp (y) = 0 otherwise. The 
right-hand side equals the Hilbert operator applied to g. The Hilbert operator is 
unitary, which yields 

/J exp (- 2tt a\J b 2 — x 2 ^ (p ( x ) | dx = J |^| 2 dx< J \<p\ 2 dx. ► 
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Radon Transform 



2.1 Properties 

Let E be a Euclidean space with the interior product (x,y) i— ► (x,y). Take a 
hyperplane H C E, choose a unit orthogonal vector u to H and denote by p the 
distance from the origin to H in the direction a>, i.e., (ou,x) = p is the equation 
of H — H(p,uj). At the same time — (lj,x) = —p is another equation of the same 
hyperplane: H(—p , —u) = H(p,u). Thus we have two-fold covering S n_1 xR^ 
A n - 1 (E) where S n_1 is the unit sphere in E and A n -\ (E) is the manifold of 
all hyperplanes in E. The topological space A n - \ (E) is homeomorphic to the 
projective space of dimension n without one point. This point corresponds to the 
infinite hyperplane in the projective closure of E. 

The Euclidean structure in E generates the volume density (measure) dV as 
well as the area density d S on any hyperplane H in V. Take an arbitrary integrable 
function / in E and define the Radon transform 

Rf(H)= [ fdS, Rf(p,u>)= [ fdS. 

JH JH(uj,p ) 

This integral is well defined for any uj and for almost all p G R. The function Rf 
is even, which means that Rf (— p, —u) = Rf (p, cu). 

Proposition 2.1 For an arbitrary f G Ti(E) the equation holds , 

(/) (<ru) = F p ^ a Rf(p,uj). 

This statement is called the “slice theorem” . 



◄ For a proof we use Fubini’s theorem: 
F p ^ <T Rf(p,oj) = f exp(-jo-p) ( f fdS) dp 

J \J {(jJ,x)=p J 



exp(— \auox)f d5dp. 



The product dS'dp is equal to the Euclidean density in E if we take an Euclidean 
coordinate system of the form (pi = (u, x) , p 2 , • • • , y n )- This density is invariant 
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with respect to all orthogonal transformations of E, hence dSdp = dx and the 
right side is equal to 

J exp(-jau>x)fdx = F (/) (<rw). ► 

Corollary 2.2 The inversion of the Radon transform can be implemented by invert- 
ing of the Fourier transform: f = F* UJ ^ X (F p ^ a Rf(p,uj)). 

Backprojection OPERATOR This is the operator that transforms an even func- 
tion g = g(uj,p) defined on S”- 1 x R to the function 

R*g(x) = \ I g((u,x) ,u>)dto. 

^ J S 71 - 1 

We put here coefficient 1/2 since the contributions of the opposite points u and —lj 
are equal. In other words the function g({u>, x ) , a;) is well defined on the projective 
space P n_1 = S n -VZ 2 and R*g(x) is the integral over this space. 

Proposition 2.3 The equation holds for any continuous integrable function f , 



R*Rf(x) = [ T y \ dy. 

r(^fi) J \x-y\ y 



◄ We can assume that x = 0 since both sides commute with translations in V. 



*Rf (0) = [ Rf{Q,uj)duj 

J s™- 1 



[ [ f(y)dSdu>= [ 5 {u<y) (f) da,, 

JS™- 1 0 J S"- 1 



JS™- 1 J{uj,y)=0 JS n ~ 1 

where dS is the area element in the plane {u,y) = 0. The Euler kernels H £ = 
^A-i/p ^ ( se e s ec tion 1.6) tend to the Dirac kernel 5$ as A —> 0. Take <j> =■ (uj,y) 
and evaluate the integral 

i,., (/> 1 ^ = fM l X ( “’ ' {y) AyAbt 

where we set z = |?/| _1 y. Apply Fubini’s theorem to the inner integral: 
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We have dou/ds = (l — s 2 ) n ^ 2 1 dv where dv is the area element in the unit sphere 
S 71-2 . Therefore this integral is equal to 



[ s (1 - s 2 ) n/S 

Jo 



/2-i ic _ r(X/2)T(n/2) |gn _ 2| 



2r ((A + n) /2) 



J s™- 1 



r(A/2)r(n/2) 

2r (A) r ((A + n) /2) 



y M 1 / (y) ay. 



The left side tends to (2.1) as A — » 0 and the first factor on the right side tends 
to 1. This completes the proof since |S n-2 | = 27r^ n_1 ^/ 2 /r ((n — 1) /2). ► 

2.2 Inversion formulae 

For n even we do the substitution p = (uj,x) + q in (2.6) and take into account 
that g^ n ~ 1 \p,u) = dp~ 1 g (p,w) is an odd function in p , 



rf(x) = lim 



*1 / 

JS n ~ 1 J\q 



9 (n ^((w^) +q,u>)dq 



9 {n +q,w) - g {n 1} (- {u,x) + q,w) 



The limit exists if g( n ^(p.u) is a Lipschitz function with respect to p. Now we 
change the order of integration and write the right side as 

J y J i + q,w) (w.x) + 9 ,w)j dw . 

Make the substitution co i— ► — uo in the second integral and see that it gives the 
same quantity as the first one. Therefore we obtain 

fix) — [ F( n_1 ) (a)— (2 2) 

} 2— 2 W 2 r(n/2) J 0 q ’ ( j 

where 

F (^) = |5^ij J 9{(u,x)+q,uj) dw 

is the normalized back projection and |5 n-1 | = 2n n ^ 2 /F(n/2) is the area of the 
unit sphere. 

Theorem 2.4 If f E 5(E) and g = Rf , then 

m = fczil! f f V<P.»VP ^ ( 2 . 3 ) 

(-j) n is— li ((w,a:) -p + 0*) n 
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◄ According to Proposition 2.1 we can reconstruct the function by means of 
two Fourier transforms: 

Rf » F p ^Rf ~ F;^ x (F p ^Rf) = f. 

Write down the right-hand side 

F L^X {Fp^aRf) = J exp(j <j{x,uj)) j exp(— j<rp)</(p,w)dpd£, (2.4) 

where we have set £ = auo and g = Rf. We pass to the spherical coordinates in 
the exterior integral. Apply the equation d£ = a n ~ 1 dadu and change the order of 
integration: 




exp(]a((x,uj) — p))cr n 1 dcr 



J g(p,u)dpduj. 



The question mark ? means that the change of variables is not legitimate since the 
integration with respect to a obviously diverges. To make it correct we introduce 
in (2.4) the increasing factor exp(— sa) and write 



/ = lim J exp(jx<rw - ea) J exp(-jap)g(p,uj)dpd^ 

= J exp (— ea + j<r((x, w) — p))a n ~ 1 da J J g(p,u)dpduj. 

Now our treatment is legitimate and we need only to evaluate the interior integral 



poo 

I(z) = / exp(— za)a n ~ 1 da, z = £ — j((x,o;) — p). 

Jo 



The integrand exp(— za)a n ~ 1 da is the trace on the chain of integration R + = 
{cr;<7 > 0} of a holomorphic form which decreases in the upper half-plane C+ = 
{Re cr > 0}. By Cauchy’s theorem we can replace this chain by the chain {cf = z~ l t 1 
t > 0}: 



I(z) 



f 



exp (—t)t n 1 dt = z n r(n) 



(”~ !) ! 

(-j) n ((x,u) -p + e'i) n 



where e’ = (2n) 1 e. Therefore 



/(*) 



(”~ !) ! 

H) n 




g(p,u)dp j 

~Tt \ 7~\ — 

((x,uu) — p + e'i) n 



The interior limit is equal to the distribution ((x, uo) — p + Oz) n . ► 

Taking the real part of (2.3) we get the expansion of the Dirac function in 
plane waves: 




2.2. Inversion formulae 
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Corollary 2.5 For any n > 1, 

Jo (x) = ( - ~ 

v ’ 2j n 



^n-1 _(CJ, 



(-I)" 1 

x — (k) n (a;, x -j- (k) n _ 



Now we obtain inversion formulas in the real form: 



Corollary 2.6 For an arbitrary integrable function f that satisfies the above condi- 
tions and g = Rf we have 



/(*) = 



u 



9 {n l) ({u,x) - p,u)dp 



for even n ( the principal value of the inner integral is taken ) and 

( 1 \( n — 1)/2 p 

for odd n . 



Remark. The integrand is an even function of u?, therefore we can remove the 
coefficient 1/2, replacing the sphere S by an arbitrary hemisphere S+. 

◄ We may assume that the function / is real- valued. Then g is real- valued too 
and we can replace the kernel in (2.3) by its complex conjugate. Taking the sum, 
we get 

2 f( x ) = ^_jy! J J R f(p > w ) [(q + 0i )~ n + (- 1 )" (q - °*) _n ] d 9 dw > ( 2 - 8 ) 

where we set q = (w,x) — p. Calculate the kernel 

Qn(q) = (q + o i)~ n + (-1)" (q - 0z) _n . 



Integrating by parts n — 1 times yields 

J Qn(q)a(q)dq = lim J [(g + ei)~ n + (-1)" (q - ez) -n ] a(q)dq. 

= ^ n ]_ I™ J [(g + si) -1 + (-1)" (q - el)* 1 } a {n ~ 1 \q)dq. 
By (1.1) we have for n odd, 

lim J [(g + ei)^ 1 + (-1)" (q - e*) -1 ] a ( ' n ~ 1 \q)dq 

= [(g + (h)- 1 - (g - 0*)" 1 ] (al”- 1 )) = = -ja ( " _1) (0)- 

Apply this equation for a(q) = Rf(q,uo) and substitute this into (2.8). This gives 
(2.7). For n even we have by (1.1) 

lim J [(q + eiT 1 + (g - ei)” 1 ] a^-^gjdg = 2 [g" 1 ] (ai”- 1 )). 

This together with (2.8) implies (2.6). ► 




34 



Chapter 2. Radon Transform 



Remark. The formulae (2.7) are local , i.e., for reconstruction of the value of / at 
a point x we only need to know the values of for hyperplanes H(p,cj) 

through the point x. Whereas the formulas (2.6) are non-local , i.e., we need to 
know f for all hyperplanes. 

Now we write both formulas (2.6) and (2.7) in a uniform way by means of 
the Hilbert operator HI. 

For any n the following reconstruction formula holds, 



f = R 



* 



H__d 

r dp 



n — 1 

Rf. 



(2.9) 



This formula coincides with (2.6) and (2.7) if we take into account that the 
operators HI and d/dp commute, M q = (— 1 ) q / 2 I for any even q and the function 
(— M) n_1 R^-^f is even in S 71 " 1 x R. 



2.3 Alternative formulae 

For even n we can rewrite (2.6) as 



(27rr/w = -(-ir /2 / srii / 



9 {n 1) (p,w)d p 



dcu 



(uj, x) - p 

= (_i)«/2 J 9 {n - 1] {{u,x) +g,o;)dg d ^ 

where the substitution p — x) + q is applied. We have 



/ 



g( n + q,u)dq 



lim 

£—►0 



/ 

J a<—i 



lim 



/ 

J q>e 



+ 

q>£ 



9 {n V((uj,x) + q,u) - g( n 1} (- (w, x) + q, -u) 



d q, 



since — q,uj) = — g( n_1 )(— { w i x ) + <7, — w) because g ( ' n " 1 \p, to) is an 

odd function. The limit exists if g ( ' n ~ 1 \p,u>) is a Lipschitz function with respect 
to p. Now we change the order of integration and we get 



(2n) n f(x) 



7 ( " +q,u) -g {n 1 \-(u,x)+q,-uj) 



d LU. 



Make the substitution uo i— > —uj in the second integral and see tha/t it gives the 
same quantity as the first one. Therefore we obtain 



\*n) Jq q 

where F(q) = - J g((u,x) + q,u)dou is the backprojection of g . 



(2.10) 
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For odd n we can commute the backproiection operator and derivatives 
in (2.9): 



Proposition 2.7 A reconstruction can be done by 

( n — 1)/2 






R*Rf, 



where the power of —A is the pseudodifferential operator 



(£)'■- 



g = F* (|£| 2p F(p) (0) 



( 2 . 11 ) 



◄ By Proposition 2.3 we have R* Rf = K*f, where K ( x ) = \ \x\ 1 . The- 

r V~ 2~ j 

orem 1.18 yields F (R*Rf) = F (K) F (/). We have F (K) = |£| 1_n by Proposition 
1.6. This together with (2.11) implies Proposition 2.7. ► 

Specify the above formulae for small dimensions. 

2D CASE. (2.9) coincides with Radon’s original formula 

(2 , 2) 

n Jo Q 

Alternative forms of this equation are: 



f(x,y) 



g' ( p , <f>) dp 



1 r r°° 

2 tt 2 J 0 ^/-oc cos<f>x + sin<f>y-p 

i f n r [a (q + p,<t>) + 9(v - p, 

2tt 2 Jo p Jo p 2 



(2.13) 



2 9 (q, 4>)] dp 



where q = cos (j) x + sin </> y. 
3D CASE. From (2.7) we find 



fix) = f s _ x 5 "((w,x) ,w)dw = --^Af s ^((w,x),w)dw. (2.14) 

Replacing integration by summation yields another reconstruction method: 



Theorem 2.8 Let f be an arbitrary function in R n ,n > 1 supported by the cube 
[0, l] n , such that D l f G Z /2 (H& n ) for |i| < n/2 + 1. If f fdx = 0 it can be recon- 
structed by means of the infinite series 

where G C Z n is the set of vectors (&i, . . . , k n ) such that gcd (fci, . . . , k n ) — 1 and 

ki > 0 . 
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Remark. Only a countable number of directions u = \k\ 1 A; is used for this recon- 
struction called “discrete” . 



◄ Prom the smoothness assumptions follows that the Fourier series absolutely 
converges: 

f( X )= E fi eX P 0 <*, x)) , /o = o. 

iez n 

Write i = qk where k E Zo and q E Z. The number q is uniquely defined for i ^ 0 
which yields 



/ (*) = E E fiexp(jq(k,x}) = EE/ / (y) exp (j<? (k, x - y)) dy 



/c£(G Q 



VJL 



fceG q 



dy 



fceG 



f(y) ,,, v 

(k,y—x)=s d{k,y } 

k { k , x) + s 



E ex P Htf s ) ds 






The function h ( s ) = Rf ^|fc| 1 k , |A;| 1 ((&, x) + s)^ fulfils the hypothesis of Pois- 
son’s summation formula (Section 1.2). Therefore we can replace the g-sum by the 
sum of Dirac functions S (s — q) taken over all integers q: 






k (fc, x) + q 

]*i’ 1*1 



► 



Remark. Glue the square [0, l] 2 to 2- torus T. The sum 



J 2 R f 



k (k,x)+q\ 

W \k\ ) 



can be thought of as the integral over the line in T that is the image of the 
straight line {y : (k,x — y) + q = 0} in R 2 . This is a geodesic of length |fc| in the 
Euclidean metric in the torus. Therefore each term in (2.15) is the average of the 
function / over a closed geodesic and the right side is the sum over all such means. 



In the general case the discrete reconstruction formula of [96] looks as follows: 



Theorem 2.9 Under the same assumption on f the equation holds 



f (*) = fo + E 



1 rDff fc (A;,#) + q 

1*1^ Vl*l’ 1*1 



fo 



where fo = J fdx. 

A proof can be done along the same lines. 




2.4. Range conditions 
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2.4 Range conditions 

Let S n_1 be the unit sphere in Euclidean space E n . It is a smooth manifold: take 
a hyperplane H in V and consider the orthogonal projection p H : S n_1 — ► H. The 
mapping pn is a local coordinate system for any point s E S n_1 where the tangent 
plane is not orthogonal to H. If </> is a function on the sphere we can extend it 
to the function <3>(x) = (j){u) that is defined in V\{0} and is constant in any ray 
x = ru, r > 0. Apply the Laplace operator A to it and take restriction to the 
sphere: 

A s </> = A3>|S n_1 , 

where A 5 is the spherical Laplace operator. A function h on the sphere is a Ca- 
binet ion if and only if the function A J s h is well defined and continuous for any j. 

Denote by 5(S n_1 x R) the Schwartz space of functions in the product 
S n_1 x R, i.e., the space of C°°-functions h = h(p,uj) in S n-1 x R such that 
\p k (p, uj) | is bounded for any i, j, k. It is easy to check that Rf E S (S n_1 x R) 

if / E *S'(R n ). We show now that the range of the Radon operator R : S'(R n ) — ► 
S(S n ~ 1 x R) is far from filling the whole space 5(S n_1 x R). First we note that 
by Fubini’s theorem 

J Rf(p, w)d p = J f(x)dx, 

hence this integral does not depend on uj\ This equation contains a continuum 
consistency condition. The image of the Radon operator satisfies many more con- 
ditions: 

Proposition 2.10 If f E S(R n ), then for arbitrary integer k >0 the function 

m k (w) = J p k Rf(p,uj)dp (2.16) 

is equal to the restriction to S n_1 of a homogeneous polynomial of degree k. 

◄ We have by Fubini’s theorem 

P k ( [ /(ajJdS') d p = [ (( u>,x)) k f(x)dx . 

\Jh(p,uj) J J 

The function £ 1 — > {^x) k is a homogeneous polynomial of degree £, hence the 
integral of the function (£, x) k f(x) is a polynomial of degree k too. The restriction 
of this polynomial to a unit sphere coincides with (2.16). ► 

This is, in fact, a complete set of consistency conditions: 

Theorem 2.11 Let g E S( S n_1 x R) be any function such that for any k = 0, 1, 2, . . . 
the moment 

p k g(p,uj)dp 

is equal to the restriction of a homogeneous polynomial of degree k. Then there 
exists a function f E S(W n ) such that Rf = g. 




f p k Rf(p,u)dp = J 
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◄ Consider the integral 

<t>(0 = J exp(— j|£| p)g (p, j|j) dp. 

The function / = F*</> will solve the equation Rf = g if we show that / £ 5(V). 
For this we need to show that </> is in the Schwartz space 5(V'). The crucial point 
is to check that 4> is smooth at the origin. Choose a number k and write 

exp(-j|£|p) = 1 + j|£|p + ^ (|£|p) 2 + • • • + h (I£|p) fe + r k (If |p) , 

where rk is a smooth function that vanishes at the origin with its derivatives up 
to the order k. Substitute in the above integral: 

4>(£) = X^l f ptg ( p > jfi) dp + / rk ^ p ^ 9 ( p > ]|[) dp • ( 2 - 17 ) 

The ith term in the sum is equal to Q|£| 2 rai(|£| _1 £). This is a homogeneous polyno- 
mial of degree i of the variables £. Therefore the first term in (2.17) is a polynomial 
of order < k. We can take any ^-derivatives of the second term up to the order 
k. Therefore the sum is, at least, in C k (V*). The function (j) belongs to C°°(V*) 
since the number k is arbitrary. It is easy to check that the all derivatives of </> are 
fast decreasing, hence 0 is in the Schwartz space. ► 

Remark. It follows that the moment conditions imply that the function g is even. 

2.5 Frequency analysis 

The range condition of the Radon transform is related to the property of its 
harmonic decomposition. The Radon transform Rf of a function / with compact 
support in a Euclidean space E n is an even function on the manifold A = lx S n_1 . 
Its harmonic decomposition is the combination of the Fourier decomposition on 
the first factor and of the expansion into spherical functions on S n “ -1 . 

2D CASE. We show that the Fourier density g decreases rapidly out of a cone H 
in the dual group E that is equal toMxZ for the parallel beam geometry. This 
means that restriction of g to H is informative enough for a reconstruction of /. 

First we find an upper bound for Fourier coefficients of the plane Radon 
transform in R 2 : 

<X<r,m) = -T J d < pj exp(— j(crs + rrvp))g(s, <p)ds. 

Theorem 2.12 For any bounded function f in E 2 with support in the unit disc B 
and any positive number e, the Fourier coefficients satisfy the inequality 

|p(cr,m)| < T exp(27r (sinh(er) |<x| - e|m|)) J \f\dV. 



( 2 . 18 ) 
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◄ We have 



g(a, m) = III, exp(— ](crp + mp))fdqdpdp , 
where p = cos p x + sin p y, q = — sin p x + cos p y. Changing variables we get 

g(a,m) = f f(x,y) [ exp(-j${x,y,<p))d(p dxdy, (2.19) 

JB LJ-ir 

where <l> = a (cos p x + sin p y) + rrup. The phase function 3> admits a holomor- 
phic continuation in (p in the complex plane C, whereas exp(— j$) has analytic 
continuation in the cylinder C/27rZ. Denote by £ = p + (mod 27 t) the complex 
coordinate in the cylinder. We have Im3> = (jgsinh^ + mip. Move the chain of 
integration in (2.19) to the circle ^ = — sign(m)e for some e > 0. By Cauchy’s 
theorem we get the same quantity 

/ 7 r /*7T 

exp(-j$(:r, y\ <p)) d<p= exp(-j$(a;, y; <p + tip)) d, <p 

-TV J —TV 

where 

\exp(-j$>(x,y,p + lip))] — exp(Im &(x,y;p + lip)) = exp (aq sinh ^ + rmp) 
and Icrgsinh^ + rmjj | < — \m\ e -f a sinh e for (x,y) G B. Therefore 

/ exp(— j$(x, y\ p + ^V ; )) dp < 27rexp (— \m\ e + asinhe) . 

-TV 

By (2.19) this implies (2.18). ► 

3D case. In this case we have the decomposition 

g(s,Lu) = ^2 / g(<7,l,rn)ew(-ias)dsYi m (uj), 

im J 

where Yi m are spherical functions. We show now that the harmonic density 
g(cr,l,m) is essentially supported in the cone £ = {a > 1} where l is the number 
of the corresponding irreducible representation of the group 0(3). We may say 
that the quantities a, /, m are related like quantum numbers of electrons in Bohr’s 
atom a > l > \m\ if we restrict ourselves to an essential support of the density g. 

Theorem 2.13 For any bounded function f in E 3 whose support is contained in the 
unit ball, the Fourier coefficients g of Rf satisfy for any positive e, k the equation 



\g(a,l,m)\ < C k 
as a 2 — > oo, 1 < a 2 < (1 — e)l(l + 1). 



1(1 + 1 )-a 2 



( 2 . 20 ) 
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◄ The spherical harmonics form an orthogonal basis on the sphere. Therefore 

g(a,l,m) = / / g(s, uj)Yi m (u) exp(-tas) dsdw, 

Js 2 Jr 

where da; = sin 0d0 A d ip is the Euclidean density on the sphere. Hence 
g(a,l,m) = J J J fY lm (0,(p)exp(-ias) dSdsdw, 



where d S is the Euclidean density in the plane p = s where p = ( uo,x ). The 
product d S d s is equal to the Euclidean measure on E, hence 



g(a, l, m) = J Ii m fdx, where I lm = J Yi m (0,ip) exp {-tap) 



da;. 



Now we estimate the integral 7/ m . The spherical functions satisfy the equation 
(A + L) Yi m = 0 where L = 1(1 + 1) and 

a 1 d . B 1 d 2 

A = - -sin 0— + 



sin 0 30 30 sin 2 0 dip 2 

is the Laplace operator on the sphere in the standard coordinates 0, cp where 
u) = (sin 6 cos ip , sin 0 sin ip , cos 0 ) . Define the gradient operator in the sphere Va = 
(a' 0 , a^/ sin#) and take into account the formula A (ab) = A(a)6+2(Va, V6)+aA6. 
For the phase function $ = ap we have 



A$ = -2$, |V$| 2 = a- 2 (|x| 2 -P 2 )- 



( 2 . 21 ) 



Write the identity 



/ 



exp(-zf>) (A + L)Yi m dw 



o — ^ 

L — a 2 \x\ 2 

and take the integral by parts. The operator A is symmetrical with respect to the 
form do; and Aexp(-z$) = - (|V<I>| 2 + 2$^ exp(-z^). Integrate by parts and 
apply (2.21) 



°-f 



Ylr 



(A + L) exp (— z4>) 
L — a 2 \x\ 2 



do; — Iim + 



/ 



Yim$ exp(-iQ) do; 
L — a 2 \x\ 2 



( 2 . 22 ) 



The last integral admits the estimate O (cr/ (L — a 2 )) , since |x| < 1 and $ = 0(a). 
This proves the assertion of the theorem with k = 1. 

Then we apply similar arguments to all the integrals in the right side of 
(2.22) and get the estimate 7/ m = O ^(<j/ (L — cr 2 )) 2 ^. Repeating this procedure 
we prove (2.20) for k = 2, 3, — ► 
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2.6 Radon transform of differential forms 



Let E be an oriented Euclidean space of dimension n and a be an n — 1-differential 
form in E whose coefficients decrease at infinity sufficiently fast. The integral 

la (H) = f a 
Jh 

is defined for an arbitrary hyperplane H = H (p, cj ) in E. We define the orientation 
of H by means of the normal vector cj ; thus the orientation of H (— p, — cj ) is the 
opposite of that of H (p, cj ). We call la the affine integral transform of the form 
a. This transform can be described in terms of vector fields. Choose a system of 
coordinates xi, . . . , x n in E and take the volume form dx = dxi A • • • A dx n . For 
a vector (tangent) field s, we set a = s h dx. Then the equation 

la ( H (p, w)) = {w, R [s] (p, w)) = ^2 u>iR [si] (p, w) 

holds, where R denotes the Radon transform, s = Sid/dxi and 

a = (— 1) 2_1 S;dxi A • • • dxi • • • A dx n . 

We have I dr = 0, if r is an n — 2-form that decreases together with dr suffi- 
ciently rapid at infinity. Therefore the form a cannot be reconstructed from data 
of hyperplane integrals. However, it is possible to reconstruct its differential. 

Proposition 2.14 If the form da is integrable, it can be explicitly reconstructed from 
data of la. 



◄ Take an arbitrary half-space U C V; let H be the oriented boundary of U. 
By the Green formula 



/ a = f 

Jh Ju 



da. 



Write da = fdx for an integrable function /. Taking the derivative of the right- 
hand side in the normal direction cj yields 



d_ 

dp 



/ a = 



L 



H(p,u) 



d a = R[f](p,u). 



Thus we know the Radon transform of /. Apply an inversion formula and we 
recover / which gives d<r = fdx. ► 

In terms of the vector field s this means that the divergence 



^ dS ( 
“ dx . 



can be reconstructed. Indeed, for a = s b dx we have da = div s dx. 




Chapter 3 

The Funk Transform 



3.1 Factorable mappings 



Let X be a Riemannian manifold of dimension n with the metric tensor g and 
Y = {Y} be a family of closed subvarieties of X of dimension k, 0 < k < n. For 
a continuous function / in X that decreases sufficiently fast at infinity, we define 
the integrals 

Mf(Y)=J^fdV(Y),YeY (3.1) 

where d V(Y) is the volume element on Y induced by the metric g. We call the 
function Mf |Y an integral transform of /. For a Euclidean space X and the 
family of hyperplanes we call this operator a Radon transform. We follow this 
terminology in any situation where the geometry is symmetric with respect to a 
transitive commutative group; another example: tori. 

We call the transformation (3.1) a Funk transform, if the manifold X pos- 
sesses a noncommutative symmetry group or in the case of algebraic varieties Y. 
For the sphere X = S 2 and the family of big circles C, this transform coincides 
with the classical Minkowski-Funk transform. Alternatively we can take the pro- 
jective plane X = P 2 instead of the sphere and the family of projective lines Y, 
which are images of big circles. 

In some special geometrical situations a reconstruction formula for the trans- 
form (3.1) can be translated for another geometry by means of the following simple 
arguments: 



Definition. Let <I> : X\ — > X2 be a smooth mapping of Riemannian manifolds 
(it need not to be an isometry), Y = {Y} be a family of smooth subvarieties of 
X\. We say that is (infinitesimally) factorable with respect to this family if 
for an arbitrary Y G Y and arbitrary point x G Y the equation holds for the 
Jacobian of <F: 



dV( 2) ($(s),g(r)) 

dV(!) (x,Y) 



= j(x)J(Y), 



(3.2) 



where dV^) (x, Z) is the Riemannian volume density on a submanifold Z C Xi, i = 
1,2 at a point x. The functions j : X M and J : Y — * R depend on $ only; 
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we call these functions Jacobian factors of the mapping 4>. An application of this 
property is obvious: the problem of inversion of the integral operator M for the 
family $ (Y) = {$ (Y) , Y G Y} is reduced to that for the family Y by 

Mf (* (y)) = [ fdV {2) = J (Y) f <D* (/) j (x) dV w 

J<f>(Y) JY 

= J(Y)M[9 * (f)j] (Y ) , 

where / is a function on X 2 and (/) ( x ) = / ($ (x)). If there is an inversion 
operator I for M |Y, we define the inversion operator for M |$ (Y) a s follows: / = 
I (/) j]). The reduction can be reverted since the mapping <F _1 is also 

factorable. Moreover, the transitivity property holds: 

Proposition 3.1 If a mapping : X 2 —> X 3 is factorable for the family $ (Y) and 
a Riemannian space X 3 , then the composition : X\ — > X% is factorable for Y 
with the Jacobian factors 



h* (x) = h (<*> (x))h (x ) , J** (Y) = J * ($ (Y)) J* (Y) . 



Example 1 . Any conformal mapping $ possesses the property (3.2) for the family 
Yk of all fc-dimensional subvarieties with the Jacobian factors j = (3>* (# 2 ) /gi ) k ^ 2 , 
J = 1, where # 1,2 are the metric tensors in 

Example 2. Let D be the unit disc in the Euclidean plane. The aiutomorphism 
of D given by G (z) = 2z (l + \z\ 2 ^ is factorable for the family of circle arcs 
A C D that are orthogonal to the circle dD. The Jacobian factors are 

2 (l + r 2 ) 1/2 



jo (z) = 



i-i*r 
(i + M ! ) 



r, Jg (A) = 



where r is the radius of an arc A. The image G (A) is the chord in D that props 
up the arc A. 



Example 3. Let E be a Euclidean space and L be a projective transformation 
of the projective closure E of E. It defines the mapping L : E\L -1 (JToo) — > 
E\L(JToo) where Hoo denotes the improper hyperplane in E. This mapping is 
factorable for the variety of ^-dimensional affine subspaces A C E and arbitrary k. 
To show this property, we take the Euclidean space E = R+E with the coordinates 
xo, • • • , x n and consider the isometrical embedding e : E — > E, e (xi, . . . , x n ) = 
(1, xi , . . . , x n ). Take a linear automorphism L of E that generates L, i.e., L = pLe 
where p : E\ {xo = 0} — ► E is the central projection, p (x) = (xi/xq, • • • , x n /xo). 
Set po (x) = xo. 

For arbitrary points 6°, . . . , b k G E denote by {6°, . . . , b k } the k x (n 4- 1) 
matrix with rows 6°, . . . , b k . The positive number [6°, . . . , b k ] is defined by 



[6 °, . . . ,b k ] 2 = 

Q=i 0 <i 1 <...<i k 
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where Bi 0)iii) i fc means the minor of the matrix {6°, . . . , b k } formed by columns with 
numbers io, . . . , ik- Note that [a 0 , . . . , a k ] = kW for arbitrary points a 0 , . . . , a k G 
E, where V is the volume of the simplex (a 0 , . . . , a k ). 



Proposition 3.2 For an arbitrary affine k -plane A in E the equation 



dV(L(x),L(A)) 
dV (x, A) 



ML (e (*)))| 



- fc -i [L(a°),...,L(a fc )] 
[a°,...,a k ] 



(3.3) 



holds, where a 0 ,. . . ,a k are arbitrary points in E that span A, i.e., 



A = jx = '^2t i a l ,^2t i = 1 j . 



◄ The quotient [L (a 0 ) , . . . , L (a fe )] / [a 0 , . . . , a fc ] does not change if we replace 
a 0 , ... ,a k by another set of vectors that span A. Therefore it is sufficient to check 
(3.3) for an infinitesimal simplex (a 0 , . . . ,a k ), i.e., for the case 



aP = x, a 1 = x + eu l , po (u 1 ) — ■ ■■ — p 0 (w fc ) = 0. (3.4) 

Suppose that the mapping L preserves the n- plane e(E). Then the quantity 
[L (a 0 ) , . . . , L (a fc )] is equal to k\W where W is the volume of the simplex with 
vertices L (a 0 ) , . . ., L (a fc ) and the right side of (3.3) is equal to W/V since 
Po (L (e (x))) = 1. This proves (3.3) for any L = Lo preserving e (E). This is true 
also for arbitrary matrix L whose first line is equal to (Z,0, . . . ,0) since dividing 
by l does not change the right side of (3.3). Now we check it for any coordinate 
permutation of the form 



L j (*^) (*^7 5 X\ ) • • • j Xj — \ , Xq , Xj -\- 1 , . . . , X n ) . 

This will imply the proposition since an arbitrary operator L is a composition of 
transformations Lo, Li, . . . , L n and the right side of (3.3) is multiplicative with 
respect to composition. 

Take j = 1 for simplicity of notation. We have 



L(l, £ n ) = ( 1 , 2 / 1 ,... , 2 / n ) , yi = l/xi,2/2 =x 2 /x 1 ,...,2/n =x n /xi. 

Take a simplex of the form (3.4). Its volume is equal to e k [x.u 1 , . . . ,u k ] (up to 
the factor (A;!) -1 ). By (1.25) the volume of the simplex (L (a 0 ) (a k )) equals 

e k [L (x) , v 1 , . . . , v k ] where 



v J = x 



-2 



^0, -u{,xiu J 2 - x 2 u { , . . ,,xi u J n 



x n u 






l,...,k. 



It is easy to check that 

{L(x),v l ,...,v k } =x 7 1 



/ Vi 1 2/2 ••• Vn \ 

u\ 0 u\ ... v} n 



\u k 0 u k ... u k ) 

= xf 1 {L (x) , L (u 1 ) ,...,L ( u k ) } . 
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Therefore 



dV(L(x),L(A)) 
dV (x. A) 



[L = ^_ fc _! [L ( x ) , L (u 1 ) , . . ■ , L (u fc )] 

[x^u 1 , . . . ,u k ] 1 [xjtz 1 , . . . , u k ] 



which coincides with the right side of (3.4) since x\ — po (L ( e (x))). ► 



Example 4. Let S be the unit sphere in a Euclidean space E with the center in 
the origin, E\ {0} — > S be the central projection. Take an arbitrary sphere S in 
E\ {0}. The projection defines the mapping n : S — ► S. It is factorable for the 
family of /c-spheres FnS where F is an arbitrary k + 1-subspace of E. The volume 
relation is 



dV (x,Ffl S) |x|* +1 

dV(x 0 ,FnS) “ |(x,i-e)r F ’ Xo ~ n ^ 



(3.5) 



where £ is the center of S', rp is the radius of the sphere FnS. The Jacobian 
factors are j (x) = |(x,x — £)| -1 |x| fc+1 , J (F) = rp. If the origin is inside 5, the 
dominator (x,x — £) does not vanish on 5; otherwise it vanishes at each point 
where a ray from the origin is tangent to S. 



Example 5. We can take a hyperplane H C E\ {0} instead of the sphere S in 
Example 4. The central projection n : H — > S is again factorable: 



dV (x, F n H) [x[ fc+1 

dV(x 0 ,Fr\S) dist(FnpO)’ x ° W(X) 



(3.6) 



This follows from (3.5) if we move £ to infinity along the line orthogonal to H at x. 

Example 6. CORMACK’S CURVES. Take a positive number a and consider the 
family of curves in a Euclidean plane E, 



A a ( P , 9) : r a cos (a (</> — 6)) = \a\ < tt/2ol, 

called a-curves. This family is rotation invariant. In the case a = 1/2, 1,2 the 
curve A a is a parabola, straight line and one branch of a hyperbola, respectively. 
For a positive (3 the (3- curve is defined by the equation 

B 0 ( p , 0) : p cos ((3 ( 0 - 4>)) = r 0 , \0\ < 7t/2/3. 

For (3 = 1/2, 1,2a /?-curve is a cardioid, a circle through the origin and one branch 
of a Bernoulli lemniscate, respectively. These two types of curves are dual: fixing 
a point x in the plane, the parameters (p, 0) of the curves A a through x belong 
to a curve B a and vice versa. For a function / in E, Cormack defined the integral 
transform for a-curves 

Mf(p,9)= [ f (x) ds 
JAfaO) 

and similarly for ^-curves. He found in [11] the inversion formulae for these integral 
transforms by means of harmonic decomposition / and Mf with respect to the 
polar angle. 




3.2. Spaces of constant curvature 
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If a is an integer or inverse integer we can apply an alternative method by 
reduction to the Radon transform. Let a = 1/n; consider the conformal mapping 
w = z n defined in E\ {0}. Any curve Ai/ n coincides with the image of a straight 
line L and 

f f (w) ds = f f (z)n\z\ n ~ 1 ds. 

JA 1/n JL 

Prom this equation we know all line integrals of g (z) = f (z)n\z\ n ~ 1 and can 
recover g by inverting the plane Radon transform. In the case a = n the image of 
A^n is again a straight line L and 

J f(z)ds = J f(n\z\ n ~^ ds. 

The inversion of the a - transform can be again done by means of Radon’s formula. 
The Cormack transform for ^-curves can be inverted in a similar way by applying 
the conformal mapping -0 (z) = z~ n , n — 2,3, — 

We consider more examples of factorable mappings in Chapter 8. 

3.2 Spaces of constant curvature 

There are three types of complete simply connected Riemannian manifolds of 
constant sectional curvature: elliptic, Euclidean and hyperbolic. A Euclidean space 
has zero curvature. Any straight line is a geodesic and vice versa. An elliptic space 
of dimension n is the real projective space P n = S n /^2 with the metric inherited 
from the unit sphere S n C E where E is a Euclidean space of dimension n + 1. The 
sectional curvature of the elliptic space is equal everywhere to 1. For a subspace 
FcEof dimension two, the intersection F n S n is a big circle; its image in the 
elliptic space P n is a closed geodesic curve 7. For an arbitrary subspace F the 
manifold Y = Fr lS n /Z 2 is a projective subspace. It is a totally geodesic manifold 
of P n , since any two points of Y can be connected by a geodesic 7 C Y. Vice 
versa, any closed totally geodesic submanifold of the projective space is equal to 
a projective subspace. 

A hyperbolic space of dimension n can be constructed in a similar way. 
Choose Euclidean coordinates xq, £ 1 , . . . , x n in E and consider the hyperboloid 
Q C E given by 

xl = x\ H b x 2 n + 1, xo >0. 

This is one sheet of the quadratic hypersurface. The hypersurface Q endowed with 

the induced pseudo-Euclidean metric dcr 2 = — dsg + dx 2 H ls a hyperbolic 

space H of sectional curvature —1. For any subspace F C E of dimension 2 the 
intersection FT 1 Q is a closed geodesic curve; for an arbitrary F it is a totally 
geodesic submanifold of H. 

Alongside this, we consider the Euclidean submanifold E = {xq = 1} and 
hemisphere S+ = {x; |x| = 1 ,xq > 0} endowed with the metric ds 2 = dxg + b 
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ds^. This is a model of the elliptic space P n . The central projection ir in E\ {0} 
defines the algebraic mapping 



Q E A S+. (3.7) 

We have ir (S+) = E and this mapping is invertible. The set n (Q) is an open ball 
of radius 1 . Bv the aforesaid, the intersection of these subvarieties with a, subspace 
F i: Ic space, 

res] 




Figure 3.1: Spaces of constant curvature 



Example 7. The mappings (3.7) are factorable, namely for an arbitrary subspace 
F C E of dimension HI we have 



dV s (y,FnS + ) 
dV E {x, FnE) 
dV H (y,FnQ) 
dV E {x, F n E) 



(1 + d 2 (F)) 1/2 (l + |a;| 2 ) -(fe+1)/2 , 
(1 - d 2 (F)) 1/2 (l - |z| 2 )” (fc+1)/2 , 



(3.8) 

(3.9) 



where n (y) = x is a point in E where the volume forms are compared and d (F) = 
distE (F (1 E,0). Note that (3.8) is equivalent to (3.5). 

For a Euclidean space E and for a hyperbolic space H and for the sphere S we 
consider the integral transform (3.1) on the variety of totally geodesic manifolds of 
dimension k. In all three cases, (3.1) is called a geodesic integral transform. We call 
the geodesic integral transform in a Euclidean space flat integral transform. The 
flat integral transform coincides with the Radon transform for affine subspaces of 
dimension k = n — 1 and with the X-ray transform for k = 1. 




3.3. Inversion of the Funk transform 
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Corollary 3.3 The geodesic integral transform in hyperbolic space H n is equivalent 
to the flat integral transform in the unit ball B C E n . The flat integral transform 
in E n is reduced to the Funk transform in the elliptic space S™ . 

The inverse reduction does not hold: the Funk transform is reduced to the 
affine transform only for functions / in S+ that vanish sufficiently fast on the 
equator, more precisely, the following condition is necessary: f (y) = o (?/q + 1 ) as 
y e S+,2/0 -> 0. 



3.3 Inversion of the Funk transform 



We write an inversion formula for the integral transform Ms acting on the manifold 
of big spheres Y C S = S n of dimension n - 1. Let S* be the dual sphere; a point 
z e S* defines the polar set z 1 = {x; (x,z) = 0}. It is the projective subspace of 
S of dimension n — 1. Let g be a bounded function on S; the function Msg (^ _L ) 
is defined in S. Take the point y = (1,0) G S C E n+1 and consider the family of 
spheres {z; (z,y) = cosf} , 0 < <j> < 7r/2 in S n . The average of Msg (^- L ) over the 
sphere equals 



G(^y) = 



1 



| S n 






(z,y)= cos (f) 



M n . l9 (z x )dS. 



Theorem 3.4 If n is even, then any sufficiently smooth even function g on S n is 
reconstructed by 



9(y) = 7^Z 



(- 1 ) 



i/2 



2n-2 7r n/2p (n/2) 



L 



V 2 / o 

COS — COS (f) 

o<p 



n— 1 



G(4) 



d 4> 
sin (f) 



◄ By (3.8) G(4>,y) = sec (j)F (y, tan^) where F (r) is the spherical mean of 

/ 2 \ _n / 2 

hyperplane integrals M E / ( H ), / = ( 1 + \x\ ) g over the sphere 

{d ( H ) = r = tan . 

By (2.2) the reconstruction is given by 

f°° dr 

g(y) = c n F<- n ~ l \r)- 

JO r 



rn/2 

= C n 

Jo 

where c~ l = 2 n_2 (— 7r) n ^ 2 T (n/2) 
Theorem 3.5 If n is odd, we have 

^^(n— i)/2 



,9 

cos (b — cos (b 
d(p 



n— 1 



G{4>) -7 



dtp 



sin 0’ 



y(y) 



2(2n) r 



-)l* r 

_1 Js n ~ 1 



cos (b — cos 6 
of) 



n— 1 



Msg (0,c<;)dt<;. 
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◄ This follows from (3.8) and (2.7). 

The integrand can be simplified as follows: 



d 



COS 4> — COS (j) 

ocp 



n — 1 



d_ 

d<f> 



n — 1 



U= 0 = l — I + c n, 2 

for some integers c n? fc, k = 2, 4, . . . , n — 1. 



dcj) 



n — 3 



+ * * * -j- C n n _ i 



3.4 Radon’s inversion via Funk’s inversion 

Funk’s inversion formula. For n = 2, integrating by parts yields 

r 71- / 2 / cog 2 (j) Q 



-ng 



(y)= [ 

Jo 



sin</> dcj) 
J cos 2 (j) 



— cos (f) j G ((f)) dcj) 



-jf( 



sin <fi 



+ sin <j> G' (<f>) d(f> 






which implies 



g(y) = ~l / 

™ Jo 



* /2 d G{y,<j>) 



+ 






. , , — \t)l ry I • (3.10) 

/o sin0 7r V 2/ 

Let / be an arbitrary continuous function in E 2 such that f -- o at 

infinity. Define the function g(y) = (1 + \x\ 2 )f(x) on S + where n(y) = x . It tends 
to zero as y approaches the equator of S+; we set g = 0 on the equator and extend 
it to S 2 as an even function. By (3.8) we have, for an arbitrary big circle C, 

M s g(C)= [ gdV s = (l + d 2 (L ))^ 2 [ /dV E = (l + d 2 (L)) 1 / 2 M E /(L), (3.11) 
Jc Jl 

where L = n (C) is a straight line in E 2 and M E is the flat integral transform in 
the Euclidean space. If we know M E /, we can determine Msg and apply (3.10) 
for the point x = 0, yo — (1, 0) taking into account that the second term vanishes. 
The quantity G(y, 0) is the average of Msg(C) for big circles whose spherical 
distance from y is equal to cj). We have 1 + d 2 (L) — cos -2 cj) for L = tt ( C ) and by 
(3.11) the right side is equal to 



;/ 

Jd(L )= t 



M E /(L)d</> 



F(x, tan 4>) 



where 



27TCOS «7d ( L )=tan 0"'“'' ' ' cos</> 

F(x,r) = J-f MeJ (L) d(f>. 

Jd(L)=r 



(3.12) 




3.5. Unified form 
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Substitute this equation in (3.10) and change the variable </> to r — tan</>: 



1 f°° d(sec 4>F(y,r)) 
n Jo sin (p 



We have 

d(sec 4>F(y,r)) d F d F 

Y )J ~ = — + rdF + Fdr = — + d(rF). 
sin <p r r 

The last term vanishes after integration over the ray (0, oo) since the product 
rF vanishes at the ends. The equation rF = o(r _1 ) for r — > oo follows from 

g = o ^|x| -2 y This yields 



/ (0) = g(y) 



ir 

7T Jo 



d F (0,r) 



Moving the origin to an arbitrary point xgE, yields 

dF (x, r) 



= f 

n Jo 



This is Radon’s formula, see Section 2.2. 

The same arguments applied to the projection n : Q 
(3.9) the inversion formula for the hyperbolic plane 




dG (x,q) 
sinh q ’ 



(3.13) 



E give by means of 



(3.14) 



where G is again defined by (4.12) and q is the hyperbolic distance. 

Comparing the formulae (2. 10), (3. 13) and (3.14), we see the obvious similar- 
ity. The form of the denominators: r, sinr, sinhr shows direct impact of Euclidean, 
elliptical and hyperbolic geometries, respectively. 



3.5 Unified form 

Consider a vector space V = V n+1 with a coordinate system x = (xo, . . . , x n ); 
define the volume density dx = dxo A - • • Adx n . Take the Euler field £ = ^2 xjd/dxj 
and define the even differential form 

n 

u = e b dx = ^ (— iy Xjdxo A • • • dxj • • • A dx n . 
o 

Let / be a homogeneous function in V\ {0} of degree — n, i.e., / (tx) — t~ n f (x) 
for any xGE,x/0,t>0. Define its Funk- Radon transform (FR-transform) as 

I(f)(0 = \j K f(x)H(tx))u(x), (3.15) 

where K is a hypersurface in V\ {0} that meets each ray with the source at the 
origin only once and £ runs over the dual space V*, £ ^0. 
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Proposition 3.6 The integral (3.15) does not depend on the hypersurface K that is 
homotopic in V\ {0} to the sphere S centred at the origin. It is an even homoge- 
neous function in V*\ {0} of degree —1. 

◄ It is sufficient to check that the integrand is a closed form (current). The 
form dx is closed and d (e b dx) = (n + 1) dx which yields d (fS^uj) = d (fS%) A 
uj + f5% (n+ l)dx where S % (x) = S ((£,x)). On the other hand, e (uj) = 0 and 
e ( fS ^) = — (n + 1) /($£, since the Dirac function 6% is homogeneous of degree —1. 
This yields ehd (fS^uj) = £ (fS%) uj + (n + 1) f5% (e h dx) = 0 which proves the 
first statement. To check the second one we note that the Dirac function 5^ is even 
and homogeneous of degree —1 with respect to £. ► 

Take an arbitrary even function g on the sphere S and extend it to a function 
/ in V\ {0} as a homogeneous function of degree — n. Then 

I [f](t) = l [ 9(y)ut= t g(y)u(, 

z J sn^-L Js+nz-L 

where is the polar set {y : (£, y) = 0} and a = uj/d(£,y). Introduce the Eu- 
clidean structure in V by means of coordinates x, then the restriction of uj to 
the oriented unit sphere S coincides with the area density dx/dr. If £ is a unit 
covector, uj^ coincides with the Euclidean area density on This implies that 

/ coincides with the Funk transform. 

To compare (3.15) with the Radon transform, we take the Euclidean hyper- 
plane Ei = {y : 7/o = 1} in E. If h is a function on Ei, we can extend it to E as 
an even homogeneous function of degree — n as follows: 

f(y) = \yo\~ n h(y 1 /\y 0 \,...,y n /\yo\)- 

We have w\Ei = dxi A • • • A dx n and I [/] (£) = R [h] (H) where H = E\ fl 
is a hyperplane in E\\ H = H (p,0), where 6 = (£, rj)r] — £, p = (£, rf) / 10| and 

7? = (1,0,. ..,0). 

Similarly the geodesic transform on the hyperbolic space is reduced to the 
transform (3.15) for functions / supported in the cone 

[y -yo > G/i + •■• + &£) 1/2 }- 

Corollary 3.7 The equations (3.8) and (3.9) hold. 

Inversion formulae. Let uj * be the n-form in E* defined similarly to uj. Take 
an arbitrary half-space E+ in E. 

Theorem 3.8 The operator I maps the space of even homogeneous functions f in 
E of degree —n isomorphically to the space of even homogeneous functions g in 
E* of degree — 1 . The inversion is given by the transformation 

J[9](x) = —^— I [ 0(O* (n “ 1) (M)"*(O 

(27TI) Jk* de + 




3.6. Punk- Radon transform and wave fronts 
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for odd n and 



J[9} = 



i(n — 1)! f 

(2m) ri ~ 1 Jk*h e. 



0(0 M “ n ^* (0 



for even n, where K * is an arbitrary hypersurface in E*\ {0} that is homotopic to 
the unit sphere. 



3.6 Funk-Radon transform and wave fronts 

Definition. [46] Let u be a distribution in an open set A c V. The wave front 
WF ( u ) is the subset in A x Vq\ {0} defined as follows: a point (x, £) , x G A, £ e 
V* does not belong to WF(u) if £ = 0 or there exists a neighborhood U of x 
and a conic neighborhood V of £ in V* such that for any function <j> € D ( U ) the 
Fourier transform of cf)u fulfils 

F (4>u) (rj) = O (M"*) > »7 S V, to| -> oo, 9 = 0, 1, 2, . . . . 

The wave front set is a closed conic subset of A x Vg\ {0}. The wave front of a 
generalized function v is defined to be WF(v dx) where dx is a smooth volume 
form in V. 

This definition is generalized for distributions on a smooth manifold A. The 
role of the space X x V* is played by the cotangent bundle T* (X) of X (see 
Chapter 1). The wave front of a distribution u E D f (X) is a conic subset of Tq (A); 
the bottom note means that only non-zero covectors are collected. A conic set is 
a union of rays, so one can reduce the cotangent bundle to the set of rays in it. 

Definition. An oriented contact element in a manifold X is a set of covectors 
{(x,t£),t € M+}, i.e., a ray in the cotangent bundle. The manifold of all oriented 
contact elements is called the contact bundle of A; we denote it by C* (A). It 
is isomorphic to Tq (A) /R+ where M+ is the multiplicative group of positive 
numbers. For any submanifold Y C A there is defined the submanifold A* (T) C 
C* ( A), called conormal bundle to Y; it consists of all pairs (x,u) where x G Y 
and u is orthogonal to Y. 

The wave front of a distribution in a vector space V can be easily defined 
in terms of the Radon or the Funk transform. The Funk transform helps to do 
it in a more natural way. We embed the vector space V to the unit sphere S in 
a Euclidean space E of dimension dim V + 1 by means of the central projection. 
Let S* be the unit sphere in the dual Euclidean space E*. We have the duality 
pairing S x S* — > R, (x,£) i— > £x. Let F C S x S* be the hypersurface where 
(£, x) = 0. The manifold F can serve as the contact bundle C* ( S) over S : for a 
pair (x,£) the vector £ is interpreted as a unit cotangent vector at x. The same 
manifold is equivalent to the contact bundle C* (S*), so we have the involution 
L : C* (S) — > (7* (S*) , (x,£) i— > (£,x), called the Legendre transform. 
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Proposition 3.9 For a distribution u G D' (U), U C S, a point (x,£) G F does not 
belong to the set WF (u) U —WF(u) if there exists a test function (j) in S, wiht 
4> (x) ^ 0, such that the function M [cf)u] is smooth in a neighborhood of £. Vice 
versa, if the point (x, £) does not belong to this set, then there exist a neighborhood 
U of x and a neighborhood V of £ such that M [4>u] is smooth in V for any test 
function 0, supp 0 C U. 

The notation — G means the set of points ( x , — £) for (x,£) G G. 

◄ We write the Funk transform of a generalized function on the sphere in the 
form 

Mv(£)= [ vdS, £gS*, 

•V 

where £-*- is the polar set of a point £ and dS is the area form in £-*- . We can shift 
to the Radon transform by means of (3.8), which yields 

Rw (p,uj) = (1 +d 2 (C- 1 )) 172 M [<f>v\ (£) , 

/ 2 \ -n / 2 

where w = ( 1 + |x| J cj>v and the arguments (p,cu) of the Radon transform 

relate to the point £ = (£o,£0 G S* as follows: p = £o/|£ r |,^ = £ 7 /|£'|. The 
function M [</>v] is smooth in a neighborhood of a point £ if and only if the Fourier 
transform of the left side, 

F ( w ) (;n) = J exp (-jpt) Rw ( p , w) dp, 

decreases rapidly as rj = tu tends to infinity in a conic set that contains a neigh- 
borhood of the vector uj = £'/ |£'| G E. This means that (x,u>) belongs to WF [cj)v\ 
for no x and vice versa. ► 

Proposition 3.10 If u is a distribution in S with compact support, then the set 
WF (Mu) U —WF (Mu) is contained m E = L (WF (u)). 

Remark. It follows that the singular support Mu is contained in the set 
pL(WF (u)) where p : (7* (S*) — > S* is the natural projection. In other words, 
the singular support of Mu is contained in the set of points £ G S* such that 
(x, £) G WF (u) (£, x) for some x G S. 

◄ Take a point (x,£) G T,\WF(u) and show that the point (£,x) does not 
belong to WF(Mu) U -WF (Mu). According to Proposition 3.9 it is sufficient 
to find a test function -0 in S* such that 0 (£) / 0 and the function M$ [ipMu] 
is smooth in a neighborhood of x, where M$ is the Funk transform of functions 
on S*: 




3.7. Integral transform of boundary discontinuities 
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By Proposition 3.9 there exists a neighborhood U of x such that M # [ipM [</>u]\ is 
smooth in S for any test function cp supported in U and for a test function ip in 
S* such that ip (£) ^ 0. Take for (p a function that is equal to 1 in a neighborhood 
U\ of x. Then we have 

[tpM [<fm]} + M* [ipM [(1 - <f >) «]] . 

The first term in the right side is smooth everywhere. The second term is smooth 
in U\ since the singularity of the kernel [ipM]is contained in the diagonal in 
S x S and the distribution (1 — <j>) u vanishes in U\. Therefore the sum M # [ip Mu] 
is smooth in U\. ► 

Example 8. Let x be the indicator function of a domain D C S and g be a smooth 
function in S. The Punk transform of / = g\ is equal to the integral 

M/(0= f gdS. 

JDnz 1 - 

The wave front of / is contained in the cotangent bundle N * (T) of the boundary, 
moreover it coincides with this bundle if g is nowhere flat in T. 



3.7 Integral transform of boundary discontinuities 

Let D be a domain in E n with the smooth boundary T = dD. Take a smooth 
function (p in E n that is negative in D and positive in the complement and d<p ^ 0 
in T; set <p± = max (±0,0). Consider the generalized functions 

/ A (x) = a ( X ) , A g C, (3.16) 

where a G D (E n ) is a smooth function. It is supported in D and smooth ex- 
cept for the boundary dD. Moreover, the wave front of f\ is contained in the 
conormal bundle N* (dD). According to Section 1.6 the family {f\} has analytic 
continuation [f\] in the complex plane A G C, in particular, 

[f-k-i] ( P ) = sf ] ( ap ) , k = 0, 1, . . . , p e D (E) , 

where 




is a derivative of the Dirac function on the boundary. The Punk-Radon transform of 
fx is defined in the usual sense for any A, Re A > 0. It admits analytic continuation 
too. We shall show that the singularity of this transform has singularity of a similar 
form on the dual hypersurface. Take embedding e : E n — > E n+1 , x h-> (1, x) and the 
projection p : E n+1 —> S n . The composition pe : E n — > S n is the standard central 
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projection as in Section 3.2. The Radon transform in E and the Funk transform 
in S are related as follows: 



R[fx}(p,uj) = (l + d 2 (^)) 1/2 M 





They have singularities of the same form. The support of these singularities is the 
hypersurface T* = L ( N * (pe (T) D U)) C S*, called the Legendre dual hypersurface 
to T. The Legendre dual to T* coincides with T since L is an involution. 




Figure 3.2: Legendre duality in spherical geometry 

We consider now the relation of geometries of dual hypersurfaces in Euclidean 
coordinate systems: 

Proposition 3.11 Let T be a smooth hypersurface in E n such that the Gauss cur- 
vature does not vanish at a point xo £ T. Then the dual hypersurface T* is smooth 
at the point £o that is an exterior normal to D at xo and its Gauss curvature does 
not vanish. Moreover, the signatures of curvature forms of V and of F* coincide if 
we choose properly the side D * of T* . In particular, D * is convex if D is convex. 
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◄ Choose the coordinates x in E n in such a way that (£oj#o) > 0 and replace 
£o by Denote by y = (t,x) the coordinates in E n+1 and by £ = (r, f) 

the coordinates in the dual space (E n+1 )*. The vector £o = ( — l,£o) ls orthogonal 
to (1, x 0 ) and to the tangent plane T (T). Therefore it is orthogonal to the tangent 
plane T ( G ) to the cone G C E n+1 generated by T. Take the function $ (y) = 
t 2 (j) (x/ \t\) in E n+1 ; it is homogeneous of degree 2 and coincides with </> in e (E n ) . 
Consider the equation 

V$> (y) = £ (3.17) 

for a point f G (E n+1 )* close to f 0 . For £ = £ 0 it has the solution y 0 = (l,x 0 ). 
Consider the Jacobian n -f 1-matrix of this system at this point: 



vL*(yo) = 



-2 (x 0 , V4>) + t x 0 V 2 (f)Xo 
V0- V 2 0x o 



V'0 - t x 0 V 2 </> 
V 2 </> 



where V 2 </> = (xo) , i = 1, 2 and *xo means the row vector xq. 
Lemma 3.12 The matrix V 2 3> (yo) is non-singular. 



◄ Multiply the columns with numbers 2 up to n = 1 by the coordinates of 
vector xq and add to the first column. Then multiply the lower n lines by the 
same numbers and add to the first row. This results in the matrix 



J = 



0 V'</> \ 

V</> V 2 (j) ) ' 



It is sufficient to show that the matrix J is non-singular. Suppose that J (v 0 , v) = 0 
for some n -b 1-vector (uo, v). Taking the first row, we find that (V</>, v) = 0 which 
means that v belongs to the tangent space T Xo (T). Then V(j)v o + V 2 4> v = 0 and 
consequently t uV 2 4>v = 0 for an arbitrary vector u G T Xo (T). By the assumption 
the restriction of the quadratic form V 2 0 to T Xo (T) is non-singular, which implies 
that v = 0 and Vo = 0, q.e.d. By the implicit function theorem the equation (3.17) 
has a unique smooth solution y = y (f) in a neighborhood V of the point £ 0 
such that y (£ 0 ) = yo- We have y (A£) = for A > 0, since V<£ is homogeneous 
of degree 1. Taking derivatives of (3.17) yields V 2< F(y(£)) y' (£) = /, where I 
denotes the unit matrix. This implies that 

y'(0 = (V 2 $(y(e)) _1 - (3-18) 

Define *!>(£) = $(y(£))> ^ ls a ^ so homogeneous function of degree 1. We have 
vj> (£) = $ tx) = tcj) (x) = 0 if £ is the exterior normal to G at a point y G dG 
close to yo; in particular, (£o) = 0. The equation (£) =0 defines a smooth 
cone G* in V*, since V’J' (£) = V# (y (0) y' (0 7 ^ 0 by virtue of (3.18). 

Consider the embedding e* : E n * — + (E n+1 ) * by e* (£) = (—1,0 and set 
r* = G* n e* (E n *). The point 0 belongs to T*and the vector y 0 is normal to G* 
at 0, s i nce yoV'h (£ 0 ) = yoV<f> (yo) = <&(yo) = 0 and since $ is homogeneous. 
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Therefore the vector xq is normal to T* at £o- Let K* be the side of G* in (E n+1 )* 
for which yo is the outward normal. The point rj = £o — £ (0, £) belongs to K* for 
small e > 0 and we have ( 77 , (0, xo)) = 1 — e |xo | 2 > 0 for small e. Further 

^ 7 ) = $ (y (ri)) - $ (y (Co)) = -e (V# (y (7?)) ,x 0 ) + 0 (e 2 ) = -s + o (s 2 ) < 0. 

This implies that < 0 in K * and > 0 in the complement to K*. Differentiating 
the equation (£) = $ (y (£)) yields 

v 2 * (0 = y' (C)‘ v 2 $ (y (0) y' (0 + v$ (y (£)) y" (£) • 

By (3.17) we find V<I> (y (£)) y" (£) = t£ y" (£) = 0, since y' (£) is homogeneous in 
£ of zero degree. Taking in account (3.18) we get 

V 2 ^(0 = (V 2 $(y(0))* (3.19) 

where the symbol * means the dual quadratic form. Denote by V 2 </> the restriction 
of the quadratic form V 2 3> to the tangent plane of T: 

T (r) = {dxeT (E n ) ; (V</>, dx) = 0} . 

The pull back of this plane by the mapping x = x (— 1, £) coincides with the plane 
(V$(l, x' (— 1,£)) ,d£) = 0. The last plane is just the plane T(r*) tangent to 
T*, since V3> (l,x'(— 1,0) = V^(£) where ^ (0 = ^ (— 1,0- Therefore (3.19) 
shows that the form Vp*^(0 coincides with the dual to the form V 2 0(x(£)) 
whereas (£o) is the exterior normal with respect to D * = K*e* (E*). The dual 
quadratic form has the same signature. This completes the proof. ►► 

Definition. We call ip the Legendre dual to the function <p. 

Example 9. Let <p{a) = a (x) — c where a (x) = 1/2 a^XiXj is a non-singular 
quadratic form and c a constant. Calculate its Legendre dual function. We have 
$ (t, x) = a (x) — ct 2 and V<F = (—2 ct, a tj Xj ) . For any vector £ = (£ 0 , £) , £ E M n , 
we find y (£) = (— £ 0 /2c£, Yli from the equation V<F (y) = £, where {b lj } is 

the matrix inverse to {a ZJ }. By definition, we have 

*(0 = *(y(0) = ^ 

i,j= 1 

and ip (£) = 4/ (— 1, £) = 5 (x) — l/4c, where b (£) = l/26^’£i£j is the dual quadratic 
form. 



Theorem 3.13 Let T = dD be a smooth hypersurface with the non-singular curva- 
ture form at a point xq , f\ be a function of the form (3.16) where a is supported by 
a small neighborhood U of x o . We have, in a neighborhood of the dual hypersurface 

r* to r n u, 



Mf\ (0 = (2t r) 



(n— 1)/2 a(x (0) + O(ifr (£)) 

|detV^(x(0)| 1/2 



9\i (0 + S, 
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where S is a smooth function and 

(i) 



0m (0 = (-1) 9/2 M = A + (n - 1) /2, 



T(/i+l) 

i/ q is even, where 'ip is the Legendre dual to (f>; 

r-\ fe\ ( iW 2 

(n) 9n(0 = (- 1 ) 

ifp is even; 

(iii) g, (0 = (-1) (9+1)/2 



r (/* + 1) ’ 



r(-M) 



7 T 



(iv) 



[V> (£)+ + cos (£)-]> 

if p and q are odd, and p ^ 0, 1, 2, . . ., and 

(_i Y +{q ~ 1)/2 



tt 

if p,q are odd and p > 0 is an integer. 



-In |^(0I, 



Remark. Let p — > m for a natural m in the case (iii). Then T (—p) — ► =boo, whereas 
the difference 

7T 

has the finite limit as in (iv). The second term in the right side is smooth and can 
be included in the term S. 

◄ Choose a point x 0 e T and set £ 0 = V</>(x 0 ); this is an exterior conormal 
vector to T at x 0 . The normalized vector £ 0 belongs to T* and x 0 is an exterior 
normal to T* at this point. The interior normal ray to T* at £o is given by the 
equation f = £ 0 - tx 0 ,t > 0 and i/j = t + o (t). The value of the function Mf on 
this ray is equal to 



Mf(t)= [ a(x) 

JH(x 0 ,t) 



(pi ( x ) dV 

r(A + 1) (£o,dz)’ 



(3.20) 



where H (xq, t), denotes the hyperplane {{x — xq, £o) = t} and (\V is the Euclidean 
volume form in E. We can take Euclidean coordinates t = (xo — x,£ 0 ), yi, ■ ■ ■ , Un- 1 
instead of x and replace (p by the function <p (t, y) = (pa ( y)—t . By Morse’s lemma we 
can change the coordinates y to z in such a way that <po (y) = Q{z) = 1/2 ±z? 

is a quadratic form. Then we can write 

Mf ( t ) = a (x 0 ) J h (z) d z + O (ip) , 



where h G D (M n ) is a function that equals 1 in a neighborhood of the origin; S is 
an integral of the form (3.20). Now we proceed by means of Fubini’s theorem 



1 



r(A + 1) 



J (Q(*)~ 



t)_ h ( z ) d z — 



JL 



h ( z ) d z 



T (A + 1) y J q— s d Q 



(t — s)+ d s 
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where the surface {Q = s} is oriented by the form d Q. Apply Proposition 3.14 to 
the inner integral. To complete the proof we need to evaluate the convolution 



/(A) = r 

Jo 

According to (1.7) we find 



,( 71 - 3)72 



(t-*Y 



r((n-l)/2)r(A + l) 



ds. 



7(A) = 77 a * 77 ( " -1)/2 = H». 



This implies the formulae (i) and (ii). The convolution 



J = 




s (n- 3)/2 

r((n-l)/2) ln|S ' 



(t - s ) 

r(A + i) 



is a bit more complicated and we have introduced the factor e G S (M) which is 
equal to 1 in a neighborhood of the origin. This helps this integral to converge at 
infinity. We need to evaluate this integral up to a smooth term for the case when 
the number k = (n — 1) /2 is natural. The convolution J can be written in the 
form 



J = F* 



F 



s+ds 

r(A + i) 



ok— l 



r (k) 



s ( s ) In |s| ds 



By (1.5) and (1.2), 



s+ds 

r(A + i) 



=r A " 1 (^-o»r A " 1 

= j- A - 1 (0+ A_1 — exp(— A7T?)crZ A_1 ) , 



F {rjk)' e ^ ln n ds ) = ~ k ( a + k - ex p ( fc7r *) a ~ k ) * f ( ct ) • 

We can take £ such that e is compactly supported by [—1,1]. We can neglect the 
term *£ for |cr| > 1. Therefore the product equals 

-^r A_fc (o-+ A_fe_1 + exp (— 7T (A + k) i) aZ x ~ k ~ l ) 

for |<t| > 1. Apply the inverse Fourier transform and get, for fi — A + fc^0, 1,2,... 

F (a+^da) = T (-/,) (-j)" (s + 0») M , F 

= r(-M)j^(s-0i) M , 

2 J = — j ~^F* + exp (— a - !^ -1 ) 

= T (— /i) exp (—n/ii) [(s + 0 1 )^ + (s - 0z) M ] 

= —2 T (— (i) exp (—n/n) [s+ + cos (fin) s(i] . 



This implies (iii). 
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To check (iv) for fi = m we add J (fi) = T (— /i) exp (— n/n) s 171 and take 
fj, —* m for a natural m. We get 

h + J = — T (— fi) exp (—7 T/n) + cos (fin) s^_ — s m ] 



-T (— g) (g — m) exp (— 7 rgi) 



S+ - S'\ 



g — m 



+ 



cos (gn) s^_ — cos (ran) s™ 
g — m 



In I 



ml 



since T (— g) (g — m) — * res_ m T = (— l) m 1 / 



3.8 Nonlinear artifacts 

Suppose we have corrupted projection data T ( Rf ) instead of Rf where T is 
a nonlinear function in R. Suppose for simplicity that T = T (s) is linear in 
the interval (— 00 , S), smooth and bounded for s > S (we call T the truncation 
function). If we substitute the corrupted data in the reconstruction operator of 
Chapter 2, we get a distorted image / = R^KT (Rf), K = (— <9 p H/27r) n_1 , instead 
of /. Some specific artificial details can appear in /, which we call nonlinear 
artifacts. Some of the nonlinear artifacts have clear geometrical structure, at least, 
for model functions / as in the previous section. Suppose for simplicity that n = 2 
and take the function / in the form (3.16) with A = 0, that is / = gS ((jf) where 
g > 0 is a smooth function. According to Theorem 3.13, we have locally Rf = 
A0 -1 / 2 + /i O for smooth functions h and ho where h > 0 since q = 0. The truncation 
T (Rf) does not affect the term ho if the threshold S is big enough and changes 
the profile of the singular function '0 -1 / 2 “uniformly” over T*. So the result will 
have the form T (Rf) ~ 'ip + ho where -0 has a more or less large gradient close to 
r*- If the threshold of truncation is big enough, we cau, nevertheless, recognize the 
shape of V from the local corrupted reconstruction /. More complicated artifacts 
are caused by the global structure of T that generates singularities in the dual 
curve r* . The action of truncation may be more strong in neighborhoods of these 
points. The following pictures show some simple situations. 

Fig. 3.2 shows the cardioid-like curve T and the dual curve T* shown in thick 
type. (The dual curve is obtained by means of the spherical duality and the central 
projection from sphere to the plane as in Section 3.2.) According to Theorem 3.13 
the function Rf increases faster near two cusp points of T* , where the curvature 
of the curve tends to infinity. These provoke strong line artifacts at the tangent 
lines to T at the inflection points. One more artifact line is the double tangent to 
T which is caused by amplification of Rf at the self- intersect ion point of T* , see 
Fig. 3.3. 

Fig. 3.4 shows the curve T consisting of two circles. The curve shown in thick 
type is the dual curve with respect to the origin at the bottom of the picture. 
Four self-intersection points of T* correspond to four common tangents shown by 
dotted lines. These are locations of nonlinear artifacts. 
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The truncation function T is model of a nonlinear measurement process in 
a standard CT technique: each measurement can be modelled by the function 
exp (— f fdS) where the integral is taken over a thin strip from the source to the 
detector. Therefore if the integral is big, the measurement is small and hardly 
distinguishable because of the noise, which means that big values of Rf are, in 
fact, truncated. The algorithms of type “filtered backprojection” (used as a stan- 
dard software) have the same structure R$KR with a different “regularized” kernel 
K whereas the projection data Rf is discrete and noised. Nevertheless the arti- 
facts that sometimes appear in a tomogram, have strong similarity to that we are 
discussing for the continuous model. See more details in [67]. 



3.9 Pizetti formula for arbitrary signature 



Let Q = 1/2 ^a^XiXj be a non-singular quadratic form in X = R n of signature 
(p, q). Define the dual differential operator Q* ( D ) = 1/2 Yl a ijdidj where {a^} = 
{a 1 - 7 } . Choose an orientation of X and consider the family of hypersurfaces 

X ( s ) = {Q (x) = s} , sGt with the orientation defined by d Q. Introduce the 
sequence of generalized functions in X supported at the origin 



F k M 



(2tt)” /2 Q* (D)*>(0) 

|A | 1/2 r (f + fc) 



k = 0, 1, 2, ... ; lj — (j)dx . 



Proposition 3.14 For an arbitrary smooth n-form u with compact support in R n , 
we have the following asymptotic expansions: 




3.9. Pizetti formula for arbitrary signature 
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(i) for p + q odd, 



f - 

J X(s) dQ 



k> 0 



,n/2-l 



S 

k\ 



+ r (s) , 



where -b k = (-1) 9/2 F k (x), if q even and -b k = (-l) p/2 F k (x), if P 
even ; 

(ii) for p and q even, 



( - 

J x (s) dQ 






E ^ s " /2_1 

fc >0 



A;! 



+ r (s) , 



where b% -b k = (-l) q/2 F k (x); 

(iii) for p,q odd, 



f iL 

J x(s) dQ 



^ ~2b k s n/ 2 1 ln\s 
k>0 



k\ 



+ r(s), 



where bk = r~ 1 Fk (x) and r is a smooth function which vanishes 

at s — 0 together with all derivatives. These expansions can be derived term 
by term. 
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◄ For a proof we apply the stationary phase method to the integral with the 
quadratic phase function: 




e P,g ft) Y'i Q* (d) 

|£A| 1 / 2 “j!\ j* 



%( 0 ) + 



Rn (v, t ) 

£JV+(tj+1)/2 ’ 



(3.21) 



where e p?g (t) = exp (sgn£7rz (p — #) /4) and the remainder Rn can be uniformly 
estimated in terms of the 2N + 1th derivatives of <fi. By Fubini’s theorem we can 




and we can recover the Funk transform of the Leray form by applying the Fourier 
transform to (3.21). First we multiply (3.21) by 1 — h where h E D (R) is an 
arbitrary function that equals 1 near the origin; the singularity of (3.21) at the 
origin is not relevant to the singularity of the Funk transform at s = 0. Applying 
the Fourier transform yields: 



[ J!L 
Jx{s ) dQ 



(£>))>( o) 



exp (-]ts) 



(1 h (t)) e p , q ( t ) 



R(s ), 



where R E (7 7V + n / 2 1 (R). Calculating the Fourier integrals by means of Examples 
10, 11 of Section 1.4, we complete the proof. ► 




Chapter 4 

Reconstruction from Line Integrals 



4.1 Pencils of lines and John’s equation 

Let E be a Euclidean space of dimension n and / be a function in E that decreases 
at infinity in such a way that any line integral of / converges absolutely. It is 
sufficient to assume that f (x) = O ^|x| _1_£ ^ at infinity for some £ > 0. We 
denote by g the line integral of / : 

g(x,6)= [ fds, (4.1) 

JL{x,e ) 

where L(x, 8) denotes the straight line (or the ray) through the point x E E that 
is parallel to the unit vector 8. In the case n = 2, data of the line integrals defines 
the Radon transform of / as follows: Rf (p,o>) = g (x, 8) where 8 is orthogonal to 
the unit vector u and x is any point such that (a;, x) — p. In the case n > 2, the 
manifold Ai (E) of all straight lines has dimension 2n — 2 which implies that data 
of all line integrals are redundant To avoid redundancy we state the reconstruction 
problem as follows: 

Problem: to find a reconstruction formula g\Y. i— ► / for functions / supported by 
a compact set K C E for a submanifold I C Ai (E) of lines (rays) of dimension 
n. We call such a manifold a pencil . The data of g (L) for a pencil of lines have 
no dimension redundancy. The pencil Y. should, of course, fulfil the completeness 
condition for K (see Chapter 6). 

The case n = 3 is important for applications to X-ray tomography. There are 
several classes where explicit reconstruction formulae are known: 

1. Choose a plane H C E and consider the pencil X# of straight lines that are 
parallel to H. Take an arbitrary plane H ' that is parallel to H. Any line L C H' 
belongs to X#, hence we know the line transform g(L). Apply the inversion of the 
Radon transform in H' and reconstruct the function / : H f — > C for each plane 
H' that is parallel to H. 
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2. Take a curve C C P(E) such that any plane H C E has non-empty intersection 
with C at infinity. Consider the pencil 51(C) of lines L that meet C at infinity. 
There exists a reconstruction method for this pencil. Indeed, take a plane H\ let 
c G C be a point where H meets C. Any line L C H that contains the point c 
at infinity belongs 51(C). Such lines L are parallel one to another and makes a 
foliation of H. Let T C H be a line that is orthogonal to L. By Fubini’s theorem 

[ fdS= f dt f fds, 

J H i/ T J L 

where dt is the Euclidean density in T. Thus we know the Radon transform Rf ( H ) 
for any 2-plane H in E and can reconstruct the function /. This reconstruction 
formula is given by four-fold integration. The method of the next section uses 
two- fold integration. 

3. Let T be a curve in E and 51 (T) be the pencil of rays with vertices in T. A function 
/ with compact support can be reconstructed if the completeness condition is 
fulfilled. The pencil 51 (T) is characteristic. A reconstruction can be done by two- 
fold integration. 

Note that the pencil 51# as in Class 1 is equal to the pencil 51 (C) where 
C = HD P (E) is an improper line. Class 2 can be reduced to Class 3, since the curve 
C can be transported to E by a suitable projective transformation P. The mapping 
P transforms lines to lines, hence, P(T.(C)) — 5Z(P(C)). By Proposition 3.2 a 
reconstruction formula for 51 (C) can be translated to a reconstruction formula for 
51 (P ( C )) and vice versa. 

4. Let S be a surface in E and 51(5) be the pencil of rays with vertices in 5 that 
are tangent to 5. It is also characteristic. A reconstruction of a function / with 
compact support is possible under the completeness condition. Class 3 is, in a 
sense, contained in the closure of Class 4. Indeed, take the ^-neighborhood T e of 
a smooth curve T. The surface dV £ is smooth if e is small enough. The pencil 
51 (cT e ) tends to 51 (r) as e — > 0 and the reconstruction formula of Class 4 has a 
limit which gives a reconstruction for the pencil Z (T). 

John’s equations. The manifold of lines Ai (E n ) has dimension 2 n — 2, hence 
g ( L ) is a function of 2n — 2 variables in any local chart of Ai (E) whereas / is a 
function of n variables. The function g satisfies a system of differential equations 
which express redundancy of variables. To write these equations we choose the 
coordinates x \ , . . . , x n in E and take the chart F in Ai (E) that contains all straight 
lines L that are not parallel to the plane x n = 0 in E. Write x = (x',x n ) where 
x' = (#i, . . . , x n -i). Take the point (x, 0) G L and the vector (0, 1) parallel to L; 
the vectors x,6 G M n_1 are coordinates of L in the chart F. These coordinates 
parameterize the line integral 

g(x, 0) = [ fds. 

JL(x,6) 




4.1. Pencils of lines and John’s equation 
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Proposition 4.1 The function g satisfies John’s equations 

8 2 h _ d 2 h . 

dxidOj dxjdOi ’ 



(4.2) 



where h(x, 6) = || ( 6 , 1) || 1 g(x, 0) . 



◄ The line L is given by parametric equations x' = x+t6, x n = £; the Euclidean 
line density in L is equal to d s = ||(0, 1)|| dt, consequently 



Therefore 



h(x,6) = J f(x + t6,t)dt. 
d 2 h _ J°° d 2 f 



dxidOj 



dyidyj 



x + t6 , £) Mt. 



We get the same integral formula for the function d 2 h/ dxjdOi. ► 

John’s equations (4.2) also hold for the line transform to the projective space 
P n . For this we consider / as a homogeneous function on E n+1 \ {0} of degree —2. 
The line transform g ( L ) of / is well defined for all projective lines L in P 3 (see 
Remark in Section 5.5). The function g still fulfils John’s equation in any chart 
of the manifold of lines which is the Grassmann manifold Gr 2 (E n+1 ). John’s 
equations belong to neither elliptic nor hyperbolic type. 

The reconstruction problem can be thought as a boundary value problem for 
solutions g of John’s equation where the “boundary” is our pencil Z and only the 
value of g on Z is known. If Z is the boundary of a domain A in Ai, this is the 
Dirichlet problem which is, in general, hard to study since the equation is non- 
elliptic. On the other hand, there is an obstacle for local solution of the second 
order equation, since no transversal derivative of g is known on Z. The only chance 
to solve such a boundary problem is to take a pencil Z that is characteristic with 
respect to John’s operator J = d 2 / dxidOj — d 2 / dxjdOi at each of its points. This 
is in fact, the case for all the Classes 1-4. Recall that a hypersurface given by the 
equation </> ( x ) = 0 in a space V is called characteristic with respect to a linear 
differential operator a(x,D), if a (x, V</> (x)) = 0 where <r(x,£) is the principal 
symbol of a, i.e., the sum of terms of a of maximal degree m where each derivative 
Di is replaced by (fl^)" 7 - 



Proposition 4.2 Any pencil of the Class 1 , 2, 3 or 4 is characteristic with respect 
to any John’s equation. 

◄ Classes 1 and 2 are subclasses of 3 defined for the manifold of lines in P 3 . 
It is sufficient to check the statement for the Class 4, it holds for Class 3 by 
continuity. Take a surface ScE 3 and a point p e S. Let L be a tangent line to 
S at p such that the normal curvature of S does not vanish in the direction of L. 
Then the pencil Z (S') is a smooth submanifold of Ai (E) in a neighborhood of L. 
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Take a pair of parallel planes P, Q and a coordinate system (xi,X 2 .,xs) as above 
to define a local chart in Ai (E). These planes are not parallel to L and we can 
take P through the point p. There tangent space Tl (Z (S)) is of dimension 3 and 
there is freedom for the vector (#i,02)- On the other hand, if we move the point 
p £ S in the direction 6 , the variation of coordinates x\ (L ) , #2 (L) will be of order 
o(\t\). If we move p by an arbitrary vector t tangent to S, the point L n P moves 
in a direction orthogonal to 2,0). It follows that the covector (</)[, <p' 2 , 0, 0) 
is orthogonal to the tangent space Tl (Z(5)). We now need only to check that 
the principal symbol a of John’s equation written in this chart vanishes on this 
covector. We have a (771, 772, ^1,^2) = —771^2 + 772^1, which makes the statement 
obvious. ► 

Remark. By [39] the inverse statement is true: any generic 3-submanifold in Ai (E) 
which is characteristic with respect to John’s operator coincides locally with a 
pencil of Class 4. Proposition 4.2 makes plausible that an explicit solution of the 
reconstruction problem can be found for Classes 3 and 4. It gives, however, no 
tool to find such explicit formulae. In Sections 7.7-7. 9 we meet a similar situation, 
where Darboux’s equation will be the master one. There we shall use some specific 
tools of theory of partial differential equations. 



4.2 Sources at infinity 

Let E = E n be a Euclidean space of arbitrary dimension n > 1 and P(E) be 
the projective n — 1-space of improper points of E, i.e., space of lines through the 
origin. Let S be the unit sphere in E; we have the natural surjection p : S —> P(E), 
which glues together each pair of opposite points of the sphere. Teike a curve C 
in P(E) and consider the pencil Z = Z(C) of lines in E that meets at infinity C. 
We reconstruct a function / in E with compact support from line integrals (4.1) 
for lines L e Z(C). 

Theorem 4.3 Let C C P(E) be a closed curve of the class C 1 that is not homotopic 
to a point, 0 = 6(s),0<s<S be a parametrization of C such that \6 f \ = 1. The 
formula 

•^ = ~i/ dS / ^(x + qO' {s),6(s))^- (4.3) 

gives a reconstruction for a function f E C 2 (E) of compact support from the data 
of line integrals g (L ) , L £ 51(C). The principal value of the interior integral is 
meant. 

◄ Write the Fourier transform of an unknown function 

/(£) = f exp(— j (£,x))f(x)dx, 

J E 
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where dx is the Euclidean volume density. Take a unit vector 0 that is orthogonal 
to £ and write x = y + t6, where (y, 6) = 0. Integrate consecutively over t and y £ 
0 L , taking into account that (£,x) = (£,y) : 

/(£) = [ exp(-j (£,y))y(y,0)dy, (4.4) 

Je± 

where d y is the Euclidean area element in the hyperplane 6 1 - . 

First note that the fundamental group 7 Ti(P(E)) is isomorphic to Z 2 and any 
projective line L is a generator of this group. The topological condition means 
that the curve C is homotopicaly equivalent to L. Fix a point Co G C* and choose 
the smooth curve C+ C S that starts with the point Co such that the mapping 
p : C + —> C is a bijection. The end of C+ is the point — co, otherwise the curve 
C+ would be closed in S and homotopic to a point. The mapping n transforms 
any homotopy C+ ~ pt into a homotopy C ~ pt, which does not exist due 
to the condition. Thus we have C* = C + U C_ where C_ = — C+. Take the 
parameterization 0 = 0(s) of the curve C+ such that 0(0) = co,0(S) = — cq. The 
point 0(s) = —0(s — S ) runs for s E [S', 2 S] over the curve C_. Then 0 — 0(s), s £ 
[0, 2 S'] is a C 1 -parameterization of the curve C* such that 0(2 S) — 0(0). We have 
(0 f (s),0(s)) = 0 since 0 and 0 ' are unit vectors. Choose a Euclidean coordinate 
system £ = (£ 1 , . . . ,£ n ) in the dual space E* and orientate E* by means of the 
volume form v = d£i A • • • A d£ n . 

Consider the n + 1-manifold C* x E* oriented by the form v A d s. The 
restriction of the form d (0 ( s ) , £) to C* x E* is equal to (0 ( s ) , d£) + (0' ( s ) , £) d s 
and does not vanish. Take a differential form u onC*x E* that fulfils the equation 

lu A d (0 (s) , £) = v A d s. (4.5) 

Such a form is defined uniquely up to a term x A d (0 (s ) , £) for an n — 1-form x- 
Therefore the restriction of u to the submanifold 

r = {(s,0-- (O(s),0=o, (0'(s), o>o} 

is uniquely defined since the form d(0(s),£) vanishes in T. We define an ori- 
entation on T by means of the form uj and consider the mapping 7 c • F — > 
E*, 7 c(s»£) = £• The pair II = (T, 7 c) is an odd n - chain in E*. 

Lemma 4.4 The chain II is dosed. 

◄ For any s the half-space T ( s ) = {£ : (0 ( s ) ,£) = 0, (0' ( s ) ,£) > 0} is com- 
plementary to T (s + S) = {(0' (s + S) ,£) > 0} and their orientations are given 
by the same form uj. This implies that the boundaries of these half-planes have 
opposite orientation and the lemma follows. ► 

The mapping 7 c is proper and its topological degree deg 7 c is well defined. 

Lemma 4.5 deg 7 c = 1 . 
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◄ Denote this mapping by 7c- The degree of 7 c is invariant under any de- 
formation of the curve C. We have a homotopy C ~ L in this class where L 
is an arbitrary projective line in P(E). Therefore we have deg 7 c : = deg7L- The 
mapping 7 l restricted to the interior of T is a bijection onto E*\L° where L G is 
the orthogonal complement to the plane that contains L. The form u' = v/ ( 0' , £) 
fulfils the equation (4.5), which implies that l 0 defines the orientation of T. It co- 
incides with the orientation of E since (0',£) > 0 in T. This yields the equation 
deg7L = 1- ► 

Lemma 4.5 implies that 

f(x)= ( exp(j {£,x))f(£)v = ( exp(j (£,^))/(£)7* ( v ) > (4-6) 

J e Jr 

where 7* (v) is the pull back of v under the mapping 7 = 7c- Take an n — 1-form 
ip in C* xE* that satisfies 'ip A ( 0 , d£) = v and set u = — (O', £) ip A d s. We claim 
that 7* (u) = u) in T. Indeed, the form p = v — uj defined in C* x E* fulfils the 
equation 



p A d (0 (s ) , £) = (0 f , £) (v Ads + *ip Ads A ( 6 , d£)) 
= (O' , £) (v — ^ A ( 6 , d£)) A ds = 0. 



This implies that p is a multiple of the form d (6 ( s ) ,£) which vanishes in T. We 
replace 7 * (v) by u; in (4.6) and substitute (4.4): 



Jr exp (j (£,x)) /(07* (w) 



[ ds / {0' (s) , 0 exp (j (£,x — y))rj> [ 9(y,0(s))dy. (4.7) 

Jo Jr(s) Jo(s) 1 - 



For any fixed s G [0, 25] the form — (O' ,£)'ip defines the orientation of the n — 1- 
plane T (s). We can write xp = da At, where a = (O' (s) ,£), and r is the volume 
form in the n — 2-space E* ( s ) = (0 (s) , 0' ( s)) ± . We have y — x = pO+qO'+r, where 
P = {y - x,0), q = (y - x,0'), r e E(s), which yields {£,x- y) = -aq + (r),r), 
where rj is the orthogonal projection of £ to E* (s). Substituting in the inner 
integral yields 



p pOO p OO 

/ (O', £} exp(j (£, x — y))ip = — I a exp (-jag) da / exp(j ( 77 , r)) r 

Jr(s) </o J-00 

><*(»•)■ 



(27 rq — (k ) 2 



We have J (r) dy = dg in the right side of (4.7) and y = x + pO + g#' since r = 0. 
Substituting in (4.7) gives the quantity 




1 

(27 xq — (k ) 2 



g(y,Q)dq 



^f d ‘L 



R (q — (h) 



; g(x + q0’ ,0)dq. 




